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INTRODUCTION  AND  SUMMARY 

by 

B.C.  Llzidlay 

1.1  OENERAL 

This  Final  Technical  Summary  Report  covers  the  final  stage  of  the  IBD 
progra4jme  on  MPD  power  generation  with  part  support  from  the  Advanced  Research 
Projects  Agency.  The  cbjectlves  and  achievements  in  this  energy  conversion 
programme  have  been  detailed  in  previous  Technical  Summary  Reports^^ ;  a 
bibliography  of  the  main  technical  reports  and  publications  by  IRD  staff  on  MPD 
power  generation  is  given  at  the  end  of  tMs  chapter. 

1.2  MPD  CIiOSED-CYCLE  LOOP 

The  MPD  closed«loop  facility  has  demonstrated  reliable  operation  during  this 
further  four-month  period,  resulting  in  a  substantial  amount  of  data  with  cesium- 
seeding.  Much  of  the  operation  with  power  extraction  has  been  at  duct  static  gas 
temperatures  in  the  region  of  1500^K.  Measured  plasma  electrical  conductivities 
of  about  1  mho/m  are  in  excess  of  the  calculated  values  based  on  thermal  equilibrium 
ionization^*.  Severed  watts  power  have  been  available  and  there  are  reasons  to 
expect  a  considerable  Improvement  on  this  level.  For  example,  cesium  seeding 
fractions  have  been  less  than  0.1  atomic  per  cent  (well  below  the  optimum), 
electrically-conducting  layers  have  been  present  in  the  MFD  duct,  the  electrode 
configuration  has  not  been  optimum  and  the  magnetic  field  has  been  restricted  to  IT. 

The  main  unsolved  technical  problems  in  the  loop  (and  which  are  important  in 
large-scale  concepts)  are  related  to  the  electrically-conducting  layer  which  develops 
on  the  generator  duct  walls  and  to  efficient  removal  of  cesium  from  the  helium 
downstream  of  the  generator.  However,  considerable  progress  has  been  made  towards 
effective  solutions  of  both  these  problems. 

A  continuing  programme  of  research  on  the  IRD  closed-loop  facility  is  required 
to  provide  a  full  range  of  data.  It  is  believed  that  the  current  develojaaent  of 
an  extemally-heated  MPD  duct  section  will,  by  permitting  isothermal  operation,  go 
far  towards  eliminating  thermal  boundary  layers  and  associated  effects  on  plasma 
behaviour.  Flow  boundary  layer  effects  remain  significant  and  a  facility  permitting 
much  larger  cross-section  ducts  will  undoubtedly  be  required.  Ideally,  such  a 
project  wotjild  bo  initiated  in  1965*  by  which  titip  much  more  basic  data  will  be 
available. 
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3  PLASMA  PHYSICS 


The  moat  aignificant  advancea  in  the  plaama  phyaics  programme  were  achieved 
on  the  siorowave  cavity  and  atomic  beam  experlmenta*  The  microwave  tranamisalon 
experiment  yielded  a  eet  of  meaaurementa  of  electron  concentration  aa  a  function 
of  temperature  for  a  cesium-helium  mixture  but  ran  Into  further  metallurgical  and 
mechanical  difficulties  which  are  currently  being  conected.  No  further  results 
have  boen  obtained  on  the  electrical  conductivity  experiments  described  in  the 
last  report  (IRL  64--39»  bibliography  number  45)#  effort  having  boen  transferred 
to  the  other  experiments. 
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CHAPTER  2 


LOOP  OPERATION 


Ly 

R.  Brown,  l.J.  Sla''.;or,  V.  Brown,  K«  Cass,  J.  Davidson  and  0.  Johnson 


2.1  imODUCTION 

A  flow  dlagran  of  the  present  loop  is  shown  in  Fig,  2.1,  Basic  design 
details,  manufacture  and  test  of  the  individual  loop  components  have  been  desoribed 
in  an  earlier  Technical  Summary  Report  (IRD  63“l),  Further  description  of  the 
main  aspects  of  loop  development,  the  details  of  33  high  temperature  runs  with 
helium  alone  and  20  cesium«>seeded  power  generation  runs  is  given  in  Teohnical 
Summary  Report  IRD  6i«.-3^  • 

During  the  20  cesium-sseded  runs  in  March  1964  the  efficiency  of  oeslutL  removal 
from  the  heliiaa  flow  was  low  (<2  per  cent).  Cons  ^uently,  it  was  neces.  azy  to 
completely  atrip  the  closed  loop  circuit  using  oxidation  and  alcohol  dissolution 
techniques.  Approximately  one-third  (lOO  ga)  of  the  jssiw  ‘  hroughput  was  located 
in  the  precooler  outlet,  the  remainder  being  dispersed  in  the  lubricating  oil  of 
the  helium  circulator.  No  deterioration  of  any  components  due  to  the  presence  of 
cesium  was  apparent.  The  closed  loop  facility  was  re-assembled  with  modifications 
and  additions.  Throe  runs  (RHT  33»  34  fe«d  35)  were  then  undertaken, 

2.2  MODIFICATIONS  AND  RE-AR5aaraT.Y 

2,2,1  Helium  ciroulator 

The  opportiini’ty  was  taken  before  re-assembly  of  the  ciroulator  to  reduce  the 
end  float  of  the  rotor  from  0,037  in,  to  0,015  in,  to  minimise  internal  end  leakage 
and  improve  the  pressure  ratio  across  the  circulator.  The  reduction  was  achieved 
by  machine  grinding  across  one  of  the  large  end  faces,  reducing  its  overall  thickness 
to  give  the  required  rotor  float. 

Later  figures  for  circulator  perforsnmce  indicate  that  a  further  reduction  in 
rotor  float  will  bo  necessary  for  optimum  efficiency  of  the  unit.  On  re-assembly 
the  circulator  was  checked  and  found  to  be  leak-tight  (<10*®#i  l/sec)  under  static 
and  dynamic  conditions. 
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2.2.2  Oil  filtration 


The  activated  charcoal  oil  filtration  unit  was  found  to  be  contaminated  with  traces 
of  cesium  compounds.  The  filter  units  were  replaced  and  vacuum-outg?ssod  at  500^C. 

2*2*3  High  temperature  helium  heater 

^■^he  high  temperature  heater  with  element  and  inner  radiation  shields  of  fine¬ 
grained  graphite  was  the  subject  of  intensive  development  over  the  period  from  February 
1963  to  March  1964*  Major  development  problems  were  associated  with  the  electrical 
power  connection  to  the  top  (hot)  end  of  the  electrical  resistance  element.  The  final 
form  of  the  heater  components  used  successfully  during  the  power  generation  runs  in 
March  1S%  (RHT  32)  is  shown  in  Fig.  2.2.  The  high  temperature  heater  system  was 
completely  stripped  down  after  RITP  52*  No  deterioration  of  any  component  was  evident 
and  the  heater  was  re-assembled  and  vacuum-out gassed  to  a  central  core  tonperature  of 
2200°C.  Up  to  the  present  time  the  heater  performance  has  continued  to  be  satisfactory, 
approximately  100  hr  of  high  temperature  vacuum-out  gas  sing  and  60  hr  high  tempera  ttire 
operation  with  helium  flow  having  been  achieved  without  further  dissembly  of  the 
heater  system.  The  heater  is  in  a  highly  active  state  when  operable,  readily 
absorbing  gaseous  impurities |  during  pcriols  when  the  loop  is  closed  down  or  under 
modification  the  high  temperature  heater  is  arranged  to  be  pressurized,  to  a  little 
above  atmospheric,  with  high-purity  helium. 

2*2.4  Cesitctt  Injection 

To  allow  seeded  operation  of  the  loop  in  March  1964  an  extremely  simple  cesium 
injection  system  was  devised  and  operated.  A  veiriable  and  metered  supply  of  pure 
helium  displaced  a  constant  cesium  feed  from  a  storage  tank  to  the  loop  through  a 
hypodermic  needle*  A  more  robust,  reliable  microfeed  mochanical  piston  displacement 
system  has  since  been  developed  (Chapter  4)*  Preliminary  tests  were  carried  out  with 
a  facility  capable  of  containing  and  injecting  1  kgm  cesium  and  frran  this  a  small 
capacity  (50  cm^)  cesium  displacement  cyolinder  was  developed  (Fig.  4*1 )  employing  a 
cesium  reservoir  neaujr  tank  allowing  for  rapid  refilling  of  the  cylinder.  Up  to  the 
present  time  the  efficiency  of  cesium  removal  from  the  helium  flow  has  not  been  better 
than  12  per  cent.  Cumulative  contamination  of  circulator  oil  from  cesium  carry  over 
allows  a  maximum  of  20  seeded  rms  of  up  to  15  sec  duration  at  up  to  1  gn/sec  cesium 
flow  (a  total  cesium  throughput  of  about  300  gm)  to  be  carried  out  during  ary  one  high 
temperature  nm,  and  the  small-capacity  injection  system  is  adequate.  The  system  was 
chemically  cleaned,  charged  with  cesium,  calibrated  and  attached  to  loop  in  preparation 
for  BHT  33. 

A  by-pass  with  6  in.  vacuum  valves  was  arranged  between  the  heat-exchanger  outlet 
and  pre-cooler  inlet,  allowing  for  the  introduction  and  isolation  of  various  cesium 
removal  systems.  Initially  a  single  stage  cyclone  separator  and  an  electrostatic 
precipitator  unit  were  fitted  (Chapter  4). 
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2.2.5  Generator  duct 


To  ".omplete  the  re-assembly  of  the  p  and  allow  far  high  temperature  clean¬ 
up  with  helium  flow,  a  t«aporary  generator  duct  system  was  fitted  (Fig.  2,3), 
comprising  two  staijaleas  steel  (l8/8)  side  walls  plasma-spray-coated  with  alumina, 
high  density  reczystsUized  alumina  top  and  bottom  walls  with  3  tantalum  0*5  in.  x 
0,5  in.  flush  electrode  pairs  and  9  alumina  spray-coated  tantalum  pin  pairs.  The 
system  was  arranged  for  investigation  of  the  characteristics  of  the  generator  in  a 
subsonic  mode;  the  pin  pairs  were  arranged  to  extend  by  different  amounts  into  the 
gas  stream  to  indicate  plasma  properties  in  the  core  of  the  flow. 

2.3  RgP  33 

When  major  modifications  have  been  made  (in  particular  to  the  high  temperattxre 
graphite  heater)  it  is  necessary  to  operate  the  loop  with  high  temperature  gas  flow 
to  transport  and  clean  up  particulate  and  gaseous  impurities.  The  main  objective 
of  RHT  33  was  to  ensure  that  the  loop  was  again  reliably  operable  througiiout  a  high 
temperature  cycle,  ai*d  to  ensure  cleanliness  in  preparation  for  the  power  generation 
experiment  with  the  16  electrode  pair  boron  nitride  duct  system  (RHT  35). 

With  the  generator  duct  fitted  the  loop  was  completed;  mass  spectrometer  and 
pressure  rise  techniques  showed  overall  leakage  to  be  less  than  '-lO’^^il/sec.  The 
loop  was  checked  out  and  charged  to  1  ata  with  helium.  Cold  circulation  was  carried 
out  for  approximately  30  min  until  the  Oa  and  Ka  Impurity  Ir*  els  in  the  helium  were 
<5  ppo  (Fig*  2.4).  Power  was  applied  to  the  high  temperature  heater,  reaching  thermal 
equilibrium  at  1770°R  (nossle  inlet  indicated  helium  temperature)  after  8  hr.  The 
rate  of  helium  temperature  rise  over  this  period  wm  controlled  to  restrict  the  level 
of  desoxbed  hydrogen  to  <50  ppm.  Sleotriccd  resistance  measureMnts  taken  between 
electrodes,  pins  and  earth  throughout  the  temperature  rise  indicated  a  declining 
generator  duct  surface  resistance  (from  >1  Mohm  at  500®K  to  <100  ohm  at  1770®K), 

From  this  it  was  concluded  that  electrically-conducting  surface  layers  were  depositing 
on  the  generator  duct  inner  surfaces.  At  this  steady  temperature  condition  several 
cesium  seeded  runs  were  carnried  out,  injecting  between  0.1  and  0.5  atomic  percent 
cesium  for  periods  up  to  15  sec.  Although  the  cesium  injection  was  successful,  very 
low  powers  were  extracted  with  magnetic  fields  of  up  to  IT  (less  than  IV  indicated 
and  total  current  of  '-10  mA)  and  over  successive  nins  with  time  intervals  of 
aj^roximately  15  min  the  generator  performance  continued  to  deterioiate.  The 
generator  system  was  similar  to  that  used  during  the  cesium-seeded  runs  (RHP  32); 
with  the  singificant  reduction  in  generator  performance  it  was  decided  that  "tiieoretical 
interpretation  of  results  from  RHP  33  would  be  impossible  and  injection  of  cesium 
was  discontinued. 
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During  the  hes  ag  cycle  the  observed  pressures  along  the  generator  indicated 
a  restriction  to  flow  which  increased  with  increasing  temperature  to  reduce  the 
helium  mass  flow  from  9  pt^sec  at  800 to  4.8  gVsec  at  1770°K. 

After  2  hr  at  1770^K  the  interelectrode  leakage  resistancej  although  extremely 
low,  appeared  to  stid>iliaoj  it  was  concluded  that  the  gas  flow  was  clean  and  power 
was  reduced  and  the  circulation  discontinued. 

The  high  temperature  transfer  section  and  generator  case  were  removed  from 
the  loopf  using  helium  blow-off  techniques,  the  remaining  sections  of  the  loop  and 
the  high  temperature  helium  heater  being  isolated  under  helium  pressure. 

Inspection  of  the  generator  section  showed  that  electrically-conducting 
layers  had  deposited  on  the  wall  inner  surfaces.  Indicated  resistances  were  less 
than  1  ohn/cm  across  the  interelectrode  spaces.  Bowing  of  the  stainless  steel 
side  walls  had  occurred  (Fig.  2.5);  the  more  serious  effect,  which  must  have  occurred 
at  high  temperature,  would  explain  the  secondary  nozzle  effect  noticed  during  RHP  33. 

2.4  RHP  34 

A  new  high  temperature  tantalum  (0.02  in.  thick)  transfer  section  and  cesium 
vaporizer,  boron  nitride  nozzle  and  16  (WlO^Ta  90^)  strip  electrode  boron  nitride 
generator  duet  were  ass^nbled  (Section  5  and  Pigs.  5.1,  5.2  a,  5.2  b,  and  5*3) 
fitted  into  the  loop.  The  objective  of  RHP  34  was  to  detennine  whether  a  high 
temperature  thermal  cycle  with  helium  flow  to  a  temperature  level  lower  than  that 
of  RHT  33  could  be  carried  out  without  deterioration  of  the  electrical  resistivity 
of  the  duct  inner  wall. 

The  loop  was  evacuated  and  recharged  with  helium.  Cold  circulation  was 
maintained  until  gaseous  impurities  were  <5  ppm  (Fig.  2.6).  Power  was  applied  to 
the  high  temperature  heater  and  rapid  heating  of  the  helium  flow  was  achieve*,  without 
measurement  of  significant  desorbed  gaseous  impurities  (nozzle  inlet  temperature 
>1550^K  in  45  min).  During  the  t«aperature  rise  the  inters lectr ode  resistances 
and  electrode  leakage  resistances  to  earth  remained  similar  to  those  expected  for 
bomn  nitride.  The  1550^K  nozzle  inlet  temperature  was  maintained  for  a  further 
30  min  without  a  significant  change  in  duct  electrical  resistance  characteristics. 
Fig.  2.7  shows  the  axial  distribution  along  the  generator  pairs  at  1400°K  (nozzle 
inlet  temperature).  The  internal  surfaces  of  the  tantalum  iiigh  temperature 
transfer  section  viewed  axially  through  the  high  temperature  heater  outlet  window 
remained  highly  polished. 

The  power  on  the  heater  was  reduced  in  stages  to  zero  and  the  loop  circulation 
stopped.  Throughout  the  cooling  cycle  the  generator  duct  wall  resistances  recovered 
to  their  original  values.  The  loop  was  closed  down  under  a  positive  pressure  o^ 
helium. 


2.5  RHT  35 

The  objectives  and  tne  theoretical  appreciation  of  results  for  RIIT  35 
discussed  in  detail  in  Chapter  3*  The  loop  was  unchanged  from  RHT  34*  Cold 
circulation  was  carried  c  t  to  reduce  gaseous  inpurities  to  <  5  ppm  (Pig.  2.8). 
Dxiring  a  period  of  approximately  1 5  br  cesium-seeded  power  generation  experiments 
were  carried  out  at  three  steady  temperature  levels.  At  the  two  lower  levels 
(approximately  1310°K  and  1690°K  indicated  nozzle  inlet  temperature)  the  visual 
appearance  of  the  high  temperature  transfer  section  and  the  electrical  resistance 
measurements  along  the  generator  indicated  clean  conditions.  At  the  highest 
temperature  level  ('^1900°K)  the  inner  walls  of  the  tantalum  transfer  section  became 
dull  and  apparently  coated,  the  effect  being  coincidental  with  a  declining 
interelectrode  and  earth  leakage  resistance  along  the  generator.  Generator 
performance  during  seeded  operation  progressively  deteriorated  and,  as  it  was 
apparent  that  further  results  would  be  difficult  to  analyse  seeded  operation  was 
discontinued  and  the  loop  cooled  and  closed  down  under  helium  pressure. 
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CHAPTER  3 


POWER  GENERATION  EXPERIMEPi’S 


by 


I.R.McNtb,  R^Brown,  V,  Brown,  K.  Cass,  J.Davidson  and  G. Johnson 

3.1  INTRODUCTION 

A  detailed  description  of  the  IRD  MPD  generator  has  been  given'’ On 
14  July  1964  a  series  of  power  generation  experiments  with  cesium  seeding  of  the 
helium  was  undertaken.  This  series  differed  from  the  previously  reported  seeded 
runs^*^  mainly  in  the  generator  and  nozzle  configuration  and  materials.  To  reduce 
end  current  leakage  at  the  inlet  to  and  exit  from  the  generator  channel,  the  tantalum 
nozzle  previously  employed  was  replaced  by  a  boron  nitride  nozzle,  the  exit  region 
of  the  generator  was  lengthened  by  about  2  in.  and  the  interior  of  the  first  part 
of  the  diffuser  plasma  spray- coated  with  alumina.  The  generator  duct,  with 
dimensions  of  Ji-  in.  x  in.  x  ^  in.,  was  constructed  from  boron  nitride.  Sixteen 
pairs  of  narrow  (1/16  in.  x  ^  in.)  tantalum-tungsten  alloy  electrodes  were  used  in 
place  of  the  five  pairs  of  square  (-t-  in.  x  i  in.)  electrodes  used  in  the  first  seeded 
run^’^j  the  gap  between  the  electrodes  was  large  (as  with  the  square  electrodes)  to 
reduce  the  possibility  of  insulator  and  boundary  layer  electrical  breakdown.  In 
addition  to  the  electrodes,  five  pairs  of  1/16  in.  diameter  tantalum  pins  were  used 
to  monitor  open  circuit  voltages,  two  pairs  being  in  the  nozzle  (before  and  after 
the  throaty  and  three  pairs  in  the  generator  exit  region.  Pour  pressure  stations 
and  three  thermocouples  were  also  placed  in  the  charinel  wall;  Fig.  3»1  shows  the 
relative  positions  of  these  instruments. 

3.2  FLOW  PROPERTIES 

During  operation  of  the  loop,  two  of  the  channel  thermocouples  failed.  The 
temperature  recorded  (continuously)  by  the  remaining  channel  thermocouple,  together 
with  that  at  the  nozzle  inlet,  is  shown  in  Fig,  3»2,  Also  shown  are  the  pyrometer 
measured  heater  element  temperatures. 

The  pressures  observed  at  the  four  stations  in  the  channel  and  at  the  nozzle 
inlet  are  shown  in  Table  3«1  in  atmospheres  absolute;  also  shown  is  the  pressure  ratio 
across  the  nozzle  found  from  the  nozzle  inlet  and  first  channel  pressure  station 
measurements.  To  obtain  sonic  velocity  across  a  nozzle  the  ratio  of  inlet  stagnation 
pressure  to  sonic  pressure  : 
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p.  /i»1 

“T  =  (“^)  =  2,05  (for  helium) 

must  ue  exceeded.  Table  }»^  also  shows  the  corrected  helium  mass  flow  in  the  loop. 
Fi-om  the  quasi  one  dimensional  equations  of  continuity  and  state,  the  flow  velocity 
may  be  found: 

U  =  k2.k  m  Vp 

where  m  is  the  kgm/sec,  T  in  and  p  in  ata;  the  Mach  number  is  then  obtained  from 

M^'  =  2.19  10"'*  UVT. 

The  velocity  and  Mach  number  in  the  channel  obtained  by  this  procedure  are  shown  in 
Table  The  calculated  Mach  numoer  is  always  subsonic,  contrary  to  the  behaviour 

expected  from  th«  observed  pressure  ratio  across  the  nozzle. 

3.3  ELECTRICAL  PROPERTIES 

The  cesium  injection  circuit  was  constructed  to  give  a  seeding  rate  of  about 
0.2  atomic  percent:  this  rate,  which  is  about  one-tenth  optimum,  is  necessitated  by 
the  low  capacity  of  the  present  evaporator.  In  practice,  while  cesium  was  injected 
for  only  a  short  time  (maximum  15  sec)  the  observed  voltages  in  the  channel  lasted 
much  longer,  decaying  only  slowly  to  zero;  in  many  cases  measurable  voltages  were 
apparent  several  minuti  after  seeding  stopped.  In  some  cases,  after  the  main 
seedJ.ng  effects  had  decayed,  large  and  apparently  random  voltage  pulses  were  observed 
probably  indicating  that  cesium  was  trapped  in  the  injection  system.  Fig.  3*3  shows 
average  open  circuit  voltages  on  pins  1,2,19  and  20  for  the  first  few  seconds  of 
run  2;  the  average  voltages  vary  considerably.  Table  3*2  shows  the  voltages 
measured  in  the  decay  period  following  the  main  cesium  flow  (for  several  runs)  together 
with  the  nominal  cesium  injection  time  and  that  actually  observed.  In  view  of  the 
slow  decay  of  voltage  the  observed  seeding  period  is  lOt  well  defined;  however, 
considerably  dilution  of  the  seed  occurs,  generally  by  a  factor  of  about  two.  Pin 
2  consistently  shows  much  higher  voltages  than  pin  1 ,  reflecting  the  higher  flow 
velocity  at  nozzle  exit.  In  several  cases  pin  2  shows  step  variations  suggesting 
that  the  observed  voltage  is  influenced  by  the  different  loads  (and  hence  different 
currents)  being  switched  on  to  the  channel  electrodes. 

In  run  10  ,  when  no  magnetic  field  was  applied,  pin  2  registered  a  peak  of 
I*,  volts.  This  could  occur  if  a  residual  magnetic  field  existed  in  the  magnet,  if 
the  two  pins  were  at  different  temperatures,  or  if  they  had  different  cesium  coverage , 

The  open  circuit  voltage  varies  considerably  down  the  channel.  Fig.  3»h 
illustrates  this  for  runs  3,  1C®  and  11;  the  voltage  falls  appreciably  down  the 
channel,  a  particularly  sharp  drop  being  between  electrode  pairs  12  and  13.  It  is 


likely  that  the  change  in  voltage  occurs  as  the  gas  velocity  decreases;  however,  as 
gas  temperature  was  measured  at  only  one  point  in  the  channel  this  cannot  be  verified. 

At  low  open  circuit  voltages, r everse  voltages  were  occasionally  observed,  as 
illustrated  in  Fig.  3*5:  this  could  arise  if  end  current  leakage  occurred  (see 
Fig,  25  McNab  and  Cooper*)  but  the  observed  voltages  are  so  small  that  this  is  not 
certain i 

The  magnet  coil  currents  for  each  run  are  given  in  Table  3«3;  Fig,  3.*n  shows 
the  calibration  curve  for  the  magnet.  This  curve  is  obtained  in  the  absence  of  the 
nozzle  and  generator  channel  assembly  and  therefore  with  a  cold  (or  self-heated) 
magnet.  When  the  MPD  loop  is  operating  the  magnet  core  is  considerably  heated  by 
radiation  from  the  hot  regions;  no  thermocouples  were  attached  to  the  magnet  during 
this  run  so  that  the  extent  of  the  heating  is  not  known*  An  increase  in  the 
operating  temperature  of  20C°C  will  decrease  the  magnetic  field  by  about  5  percent 
(Bozorth® ) , 

To  obtain  the  electrical  conductivity  of  the  plasma  in  the  -hannel  on  seedir^, 
each  of  the  sixteen  electrode  pairs  was  switched  on  to  five  load  resistors  (18,  68, 

280,  1000,  4700  ohm)  and  open  circuit  in  rapid  sequence  and  the  resulting  voltage 
displayed  on  a  24  channel  u.v,  recorder  trace.  This  sequence  takes  about  l-y  sec 
and  can  be  repeated  continuously®  From  the  deflections  on  the  u.v,  traces  the 
voltages  and  currents  for  each  load  may  be  obtained  euid  voltage-current  characteristics 
constructed.  The  slope  of  these  curves  gives  the  conductivity  of  the  plasma  when 
the  electrode  area  and  path  length  are  known.  While  this  procedure  is  tedious  and 
subject  to  some  error  due  to  ’hash"  emd  overlapping  of  the  traces,  it  enables  a  large 
amount  of  data,  to  be  obtained  for  a  short  cesium  injection  time.  In  addition  a  foui' 
channel  u.v,  recorder  displays  the  open  circuit  voltage  on  pins  1,  2,  19  and  20. 

Table  3.4  compares  measured  open  circuit  voltages  with  the  product  UBd.  In 
all  cases  the  measured  voltage  is  less  than  OBd,  the  maximum  ratio  VoeaV^UBd  being 
0.60.  Conducting  deposits  on  the  inner  channel  surfaces  may  be  responsible,  although 
these  apparently  only  become  significant  at  higher  temperatures,  as  in  run  13«  In 
runs  6,  7  and  8  the  measured  open  circuit  voltage  exhibits  a  saturation  effect  with 
magnetic  field  similar  to  that  found  previously  (Lindley  et  al*). 

Table  3*3  shows  electrical  conductivities  calculated  from  the  volt age -current 
curves.  In  many  cases  the  voltage-current  curves  exhibit  behaviour  idiich  is  not 
interpre table;  in  these  cases  the  best  linear  fit  is  employed.  The  slope  of  the 
voltage-current  is  the  internal  resistance  (r)  of  the  plasma,  which  is  related  to  the 
conductivity  by  o  =  l/rA  =  475/r  mho/m  where  r  is  In  ohm,  L  is  the  distance  between 
the  electrodes  (l,5  in.)  and  A  the  electrode  area.  Assuming  that  the  whole  channel 
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volume  carries  current,  A  is  the  channel  vridth  (-J  in*)  multiplied  by  the  electrode 
pitch  (4-  in.)*  This  assumption  yields  the  most  pessimistic  value  of  conductivity; 
the  assumption  may  be  tested  by  examining  the  conductivities  on  each  side  of  an 
unloaded  electrode.  Table  5-3  shons  tiiat  these  conductivities  are  not  appreciably 
higher  than  the  channel  average,  indicating  that  the  whole  volume  is  not  carding 
cxirrent*  The  true  conductivities  therefore  lie  between  the  values  given  in  Table 
3*3  and  values  four  times  larger.  Fig.  3.7  compares  the  average  channel  conductivitie 
with  thermal  equilibrium  values  in  the  same  temperature  range.  A  seeding  fraction 
of  0,1  atomic  percent  was  used  in  calculating  the  equilibrium  values  and  the  two 
curves  are  for  the  appropriate  channel  inlet  pressures  (the  lowest  measured  in  the 
channel)#  The  obsexnred  conductivities  are  up  to  fifty  times  larger  than  equilibrium 
at  the  gas  temperature.  To  account  for  these  valu33  merely  by  incorrect  temperature 
measurement  would  require,  in  the  highest  caces,  errors  of  at  least  300°C;  this  is 
unlikely. 

In  several  cases  the  voltage ^current  curves  were  as  shown  in  Vig.  3*8,  In 
these  cases  the  lowest  voltage  points  are  probably  in  error  owing  to  the  difficulty 
of  measuring  small  voltages  on  the  u.v,  trace.  For  a  low  load  resistor  a  small 
uncertainty  in  voltage  measurement  causes  a  large  error  in  the  derived  current  value. 

If  the  curves  are  correct  some  form  of  current  limitation  appears  to  be  occurring. 

Well  defined  concave  voltage-current  curves  were  also  observed  in  several  cases.  If 
these  are  correct  the/  indicate  electron  temperature  elevation,  although  the  same 
comments  concerning  Ihe  measurement  of  low  voltages  apply. 

The  expression  given  by  Hurwitz  et  al^  for  the  ratio  of  electron  to  gas 
temperature  in  a  segmented-electrode  generator  with  magnetically-induced  ionization 
is 


where  K  is  the  loading  factor,  X  the  ratio  of  specific  heats,  M  the  Mach  number  and 

and  the  Hall  parameters  for  electrons  and  ions.  Neglecting  ion  slip  and 

assuming  a  low  seeding  fraction,  so  tiiat  the  electronic  collision  frequency  depends 

only  on  electron-helium  atom  collisions  (q  „  =  6.1xl0~^°m^  as  in  McNab® )  solution 

e— He 

of  the  resultant  quadratic  equation  yields 


T 

^  =  0.5  +  0.5 


1 


+  8.88x10-^ 


(1-K)2M2b2T 

P 


For  run  3b  (electrode  8)  a  concave  voltage -current  curve  was  observed.  Fitting  the 
results  to  the  above  equation  (assuming  that  the  observed  conductivities  correspond 
to  elevated  electron  temperatures  given  by  Saha’s  equation)  values  of  the  inelastic 
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collision  factor  (5)  were  found*  For  the  three  load  resistors,  4700,  1000  and 
280  ohms,  the  values  of  0  were  8»45f  19*2  and  25»5*  As  the  load  resistance  decreases 
the  Ci^rrunt  flowing  through  the  plasma  increases j  for  these  thj’ee  cases  the  ratios 
S/l  were  02,  3®88  and  3*96  ruA”^  ,  possibly  indicating  that  as  the  current  density 
increases,  increased  amounts  of  impurities  enter  the  plasma. 

In  view  of  the  uncertainty  in  the  voltage-current  carves  at  low  voltages 
no  attempt  has  been  made  to  analyse  these  results  in  greater  detail. 

3.4  DISCUbSION 

Temperatures  measured  by  thermocouples  embedded  in  che  channel  wall  are  in 
error  through  conduction  along  the  leads.  While  these  errors  can  in  principle  be 
calculated,  in  practice  this  is  difficult  unless  the  complete  environment  temperatures 
are  known;  this  was  not  the  case  for  this  run»  In  addition  gas  flow  calculations 
are  required,  tn  relate  the  wall  temperature  to  the  midstream  static  gas  temperature. 
Introduction  n'f'  thermocouples  into  the  gas  stream  has  been  dismissed  in  view  of  the 
attempts  to  achieve  supersonic  flow. 

Attempts  viere  made  to  estimate  the  generator  channel  wall  temperature  profile 
by  measuring  the  change  in  resistance  of  the  individual  electrode  strips.  This 
procedure  was  carried  out  at  five  different  temperatures  (as  recorded  by  *he  single 
operating  thennocouple  in  the  channel^  during  the  operation.  Assuming  that  the 
resistivity  of  the  electrodes  follows  the  law  p  =  a  +  bT  the  difference  between  the 
measured  *hot'  and  'cold’  resistances  gives  a  direct  measure  of  the  wall  temperature. 
Fig,  3*9  shows  this  difference  down  the  channel  for  two  temperatures.  The  general 

tendency  is  for  the  temperature  to  fall  along  the  channel,  although  the  effects  of 
conduction  along  the  electrode  leads  may  vary  down  the  channel#  Fig.  ^,10  shows 
the  average  resistance  difference  against  thermocouple  measured  temperature  for  five 
different  times.  The  accuracy  of  these  measurements  is  in  grave  doubt  following 
the  discovery  of  several  cracked  electrodes  on  dissembly  of  the  channel. 

Pressure  measurements  in  the  channel  (Table  3«1 )  indicate  a  rise  in  pressui*© 
for  the  first  three  stations  followed  by  a  fall  for  the  laiit  station.  In  addition, 
the  sharp  drop  in  open  circuit  voltage  between  electrodes  12  and  13  (Fig,  3«4) 
occurs  at  the  same  point  in  the  channel,  and, although  a  pressure  ratio  across  the 
nozzle  of  more  than  2,05  was  achieved,  supersonic  Kach  numbers  were  never  found. 

This  behaviour  could  be  consistent  with  supersonic  flow  where  shock  occurred  between 
the  third  and  fourth  pressure  station.  However,  Table  3«1  shows  that  similar 
pressure  profiles  exist  when  the  flow  is  far  from  supersonic,  for  example  with  a 
pressure  ratio  of  1.22,  and  the  existence  of  cracked  electrodes  could  account  for 
the  sharp  drop  observed  open  circuit  voltage. 
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Further  experiisonte  are  required  to  clarify  the  sit’jation.  The  general 
discrepancies  between  predicted  and  observed  flow  properties  undoubtedly  follow 
from  (aasongst  otl'or  factors)  the  operation  of  a  nozzle  with  a  pressure  ratio  more 
than  seven  times  smaller  than  the  design  '•'^aluefl  The  attainment  of  only  a  low 
pressure  ratio  across  the  nozzle  follo"'s  from  the  inefficiency  of  the  compressor 
and  the  higher-than-anticlpated  pressure  losses  in  the  loop  components.  Under 
these  circumstances  the  nozzle  design  pressure  ratio  should  be  matched  to  that 
achievable  in  the  loop  and  consequently  the  following  nozzle  is  now  (17  August 
1964)  in  place  in  the  loop 


inlet  temperature 

1450°C 

inltt  pressure 

0.73  ata 

pressure  latio 

2.40 

throat  area 

2.272  cm? 

exit  velocity 

2288  ny^sec 

exit  temperature 

930°C 

exit  Mach  number 

1.13 

The  throat  area  above  is  greater  than  that  for  the  nozzle  used  here  and  the 
nozzle  has  been  cor.structed  merely  by  removing  material  fiom  the  present  nozzle 
throat.  The  above  conditions  are  evaltsated  for  pure  helium;  the  addition  of  cesium 
will  affect  the  flow  because  of  dissociation,  but  provided  the  seed  concentration  is 
low,  this  is  not  significant. 

A  further  complicating  effect  exists  in  the  present  geometry  because  of  the 
close  proximity  of  the  nozzle  to  the  generator  channel.  Since  the  magnetic  field 
does  not  cut-off  sharply,  flow  through  the  nozzle  may  be  affected  by  eddy  currents. 

The  simplest  way  to  overccane  this  is  to  increase  the  channel  length,  restricting 
^h‘»  magnetic  field  to  the  central  region. 

Leakage  resistances  were  measured  throughout  the  run  and  values  are  shown  in 
Fig,  3,11  for  electrode  pair  10,  The  resistance  decreases  as  temperature  increases, 
as  expected;  it  is  only  at  the  highest  temperature  that  the  leakage  resistance 
becomes  comparable  with  the  highest  load  resistance.  On  dismantling  the  loop  after 
operation  a  black  conducting  deposit,  similar  to  tliat  found  in  the  run  during  March ^  , 
was  found  on  the  inside  walls  of  the  nozzle  and  channel.  This  conducting  layer, 
together  with  the  cracked  electrodes  p.eviously  mentioned,  could  be  responsible  for 
several  anomalous  effects,  including  the  low  open  circuit  voltages,  Jhe  s!,arp  open 
circuit  voltage  drop  between  electrodes  12  and  13  and  the  low  and  .fidely  varying 
voltages  observed  in  runs  13,  14  and  15*  Sxcept  in  isolated  circumstances  it  is 
impossible  to  interpret  the.ie  last  results  owing  to  the  considerable  voltage  variation 
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on  each  reaistence  value*  There  As  no  pronounced  upwjtrd  or  downward  trend  in  the 
obser’ved  electrical  conductivities  along  the  channel  (Table  3*5)*  Consequently, 
remembering  the  relativsly  low  flow  velocity  and  anomalous  flow  behaviour,  it  ia 
unlikely  that  frozen  flow  occurred  during  the  nozzle  exi'ansion.  Any  higher-than- 
equilibrium  conductivities  therefore  indicate  magnetically  induced  extra- thermal 
ionization* 

The  maximmii  electrode  current  observed  during  the  run  was  about  1 5  mA , 
although  extrapolated  short  circuit  currents  ore  higher.  For  the  electrode  area 
(2  X  10"*’  this  yields  a  current  density  ol’  750  A/'m^ .  This  is  far  higher  than 
the  thermionic  current  density  obtainable  from  pure  tantalum  but  less  tiian  that 
from  pure  cesium;  Fig.  shows  values  calculated  from  the  Richardson-Dushman 

equation  and  three  experimental  values.  Unless  arc  spots  are  occurring,  the 
electxx)des  must  have  a  partial  cesium  coverage  which,  if  power  is  to  be  provided 
continuously,  must  be  continuously  ?replenished.  Different  coverage  factors  could 
account  for  the  variation  in  power  output  and  coriductivi ties  along  the  channel, 

3*5  CONCLUSIONS 

Many  complicated  and  interrelated  factors,  including  flew,  materia]  and 
electrical  effects,  are  at  present  combining  to  make  interpretation  of  results 
from  seeded  runs  a  complex  problem.  Nevertheless,  accounting  all  probable  errors, 
electrical  conductivity  values  higher  than  for  thermal  equilibrixim  ionization  appear 
to  have  been  observed  in  a  segmented-electrode  generator.  Further  experiments  to 
clarify  these  observations  are  planned. 
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LOOP  PARAlffiTERS  :  SEEDED  RUNS  14-  JULY 


CKSIUlt  INJECTION  PERIOD  AND  OBSERVED  VOLTAOES  ON  ELECTRODE  10 
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PRESSURE  ,THERMCXOUPLE  AND  ELECTRODE 
AND  CHANNEL 


STATIONS  IN  THE  NOZZLE 

FIG  3.1 


OBSERVED  TEMPERATURE  DURING  SEEDED  RUN 
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FIG  3.3 


VARIATION  OF  PIN  VOLTAGES  WITH  TIME 
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OPEN  CIRCUIT  CHANNEL  VOLTAGES 
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REVERSE  VOLTAGES  IN  CHANNEL 


MAGNET  CALIBRATION  CURRENT 


FIG  3.6 
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OBSERVED  VOLTAGE-CURRENT  CURVE 


ELECTRODE  RESISTANCE  DIFFERENCES  ALONG  CHANNEL 
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AVERAGt  ELECTRODE  PESISTANCC  DIFFERENCE  COMPARED  WITH 
TEMPERATURE  MEASURLC  BY  THERMOCOUPLE 
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CHAPTER  4 


CESIUH  INJECTION  AND  RECOVERY 


by 

R, Brown,  TJlrchbold,  and  A. Sands 
4.1  OPERATING  PROBLEMS 

The  ultimate  objective  for  t^-e  MPD  n-^ed-lnn-n  Is  to  continuously  supply 
and  meter  cesium  vapour  to  the  high  temperature  heliu'’  flow,  and  to  attain  the 
same  temperature  as  the  helium  by  superheating  the  cesium  vapour  before  injection 
or  through  efficient  mixing  of  the  helium-cesium  after  injection.  To  date, 
operation  of  the  loop  with  cesium  seed  has  clearly  dem  'trated  the  inefficiency 
of  the  present  vapo-dzing  and  recovery  techniques. 

It  has  been  established  that  during  injection  a  cloud  of  cesium  droplets 
form  in  the  diffuser-heat  exchanger  region.  In  recent  operation  the  cloud  density 
and  duration  has  been  monitored  at  the  heat  exchanger  outlet  station,  where  two 
opposite  windows  are  arranged  to  accommodate  a  light  source  and  photo  cell 
respectively.  The  responses  from  the  cell  are  continuously  displayed  on  a  u.v. 
recorder.  The  windows  systems  are  1  in.  ID  x  6  in.  long  and  *-  ve  remained 
completely  clean  throughout  a  series  of  runs,  including  a  throughput  of  approximately 
600  gm  cesium.  Although  the  cloud  monitoring  system  does  not  give  an  absolute 
measure  of  the  cesium  throughput  it  establishes  the  delay  interval  from  the  injection 
point  to  the  heat  exchanger  outlet  and  gives  some  indication  of  dilution  compared 
with  the  nominal  rate  of  injection. 

All  seeded  runs  to  date  have  been  associated  with  low  removal  efficiencies 
and  carry-through  of  cesium  droplets  to  the  helium  circulator.  The  original 
concept  of  cesium  removal  involved  cooling  the  helium-cesium  mixture  through  the 
heat  exchanger- allowing  the  flow  to  impinge  at  a  lew  velocity  on  a  cooled  sui*face 
at  the  base  of  the  heat  exchanger  and  finally  passing  the  mixture  through  a  cooler 
with  large  water-cooled  surface  area.  The  method  proved  totally  inadequate,  with  a 
net  removal  efficiency  of  less  than  2  per  cent.  In  the  circulator,  cesium  forma  a 
highly  viscous  complex  in  the  lubricating  oil  ,  with  subsequent  overheating  of  the 
circulator.  Resin  impregnation  was  employed  originally  on  the  circulator  to  achieve 
a  leedc-tight  unit,  so  that  the  possible  ingress  of  cooling  water  through  breakdown 
of  the  tpoxy  resin  represents  a  haiiard  and  limits  the  tolerable  overheat  of  the 
circulator.  Under  these  conditions  of  cesium  throughput  to  the  circulator  a 
maximum  of  20  seeded  operations,  each  of  15  sec  duration  and  up  to  0.5  /o  cesium 
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(total  approximately  300  gm) ,  can  be  carried  out  before  experiments  are  suspended 
to  allow  purging  of  the  circulator.  Purging  with  lubricating  oil  takes  several 
hours  of  cold  running  with  helium  alone. 

The  positive  displacement  cesium  injection  systems  employed  up  to  the  present 
for  seeded  runs  are  not  entirely  satisl'Bctory  for  short  injection  periods.  The 
metered  cesium  liquid  is  brought  via  the  injection  system  (Pig,  ^+-.1 )  stop  valve 
through  the  double  skin  of  the  loop  and  into  the  internal  vaporizer  (Fig,  4.2), 

The  liquid  path  is  relatively  long  and  on  short  term  injection  (<  1 5  sec),  the 
seeding  rate  is  not  uniform.  Continuous  open  circuit  voltage  measurements  taken 
during  seeded  runs  exhibit  pulsing  effects  and  indicate  dilution  of  the  actual 
injection  rate,  probably  caused  by  limited  and  intermittent  boiling  of  the  cesium 
vapour. 

The  problems  to  be  resolved  are  therefore  : 

to  improve  the  efficiency  of  seed  removal  and  recovery  from  the 
present  12  to  99  percent  and  so  facilitate  continuous  and  steady 
conditions  of  seed  injection  and  throughput; 

to  improve  the  efficiency  and  handling  capacity  of  the  cesium 
vaporizer  and  superheat  system  to  allow  for  seed  fractions  up  to  1 *5 
atomic  percent;  and 

to  determine  the  temperature  of  the  cesium  vapour  (probably  using 
spectroscopic  lino-reversal  techniques). 

4.2  CYCLONE  SEPARATION 

4,2.1  Introduction 

After  the  seeded  run  in  March  1964  (RHT  32)  a  cyclone  separator  was  considered 
for  the  removal  of  cesium,  the  design  conditions  being  as  follows; 

helium  mass  flow,  &  =  6,7  gn/sec 

heat  exchanger  outlet  temperature  =  100°C 

Mach  number  at  generator  exhaust  =  1 

The  cyclone  inlet  pressure  is  determined  by  the  conditions  in  the  generator 
duct  and  by  the  subsequent  expansion  in  the  diffuser.  Assuming  that  the  kinetic 
energy  of  the  gas  at  the  heat  exchanger  inlet  iz  low  compared  witn  the  kinetic 
energy  at  the  diffuser  inlet,  then 

^  ,  Ml  ^  ^  (pankhurst  and  Holder^)  (4.1 ) 

ro  2 
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Now  Y  =  1.63  for  helium 

Ml  =  1 

f  =  O465 

so  that  =  1  .82 

Po 

/rom  McNab  an'’  Cooper ,  when  B  =  C,  the  generator  duct  exhaust  static  pressure 
will  be  16000  newtoii/m^.  Thus,  the  diffuser  exhaust  static  pressure  will  be 

16000  X  1.82 


Pi  =  0.291  ata 

The  pressure  change  through  the  shell  side  of  the  heat  exchanger  wild  be 
negligible  (Appendix)  so  that  the  above  pressure  will  exist  at  the  inlet  to  t^e 
cyclone,  that  is 


Pj  =  0.291  ata 

Although  it  is  known  that  the  cesium  condenses  as  a  mist,  the  droplet  size 
is  not  known.  According  to  Perry* ,  condensed  mists  and  fogs  range  in  droplet  size 
from  0,1^  to  20^l}  it  is  therefore  proposed  to  design  a  cyclone  having  the  highest 
possible  efficiency  for  small  particles.  Since  both  pressure  loss  and  collecting 
efficiency  increase  as  the  size  of  a  cyclone  is  reduced,  the  cyclone  will  be  made 
as  small  as  possible  provided  the  cyclone  pressure  loss  does  not  exceed  0.5  psi. 

The  cyclone  will  then  have  the  maximum  efficiency.  Des^^'^  data  are  presented  in 
Perry*  for  cyclones  having  the  proportions  given  in  Fig.  4.6. 


Thus,  the  pressure  difference  between  the  inlet  and  outlet  of  the  cyclone 
is  given  as 


For  cyclones  of  the  given  proportions 
F23  =  8h2  (Peny*) 


so  that 
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AP23  =  7.405 

The  inlet  velocity  head  h2  = 

where  V2  = 

Now  &  = 

Aj  = 

= 

and,  hence  Vj  = 

and  ha  = 

from  which  it  follows  that 

dPj  3 
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As  previously  stated 

AP23  =  0.5  psi 
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and  -  9.48  cm  (=  3.73  in,) 


Using  the  proportions  given  in  Fig.  4.6  as  a  guide,  a  cyclone  was  constructed 
as  detailed  in  Fig,  4.7. 

4.2.2  Construction 

The  cyclone  must  be  inserted  into  a  six-inch  bore  pipe  downstreeun  of  the 
heat  exchanger  shell  side.  A  conical  section  changer  is  needed  to  lead  the 
inflowing  gases  smoothly  into  the  tangential  entry  duct.  This  is  shown  as  an 
integral  part  of  the  assembly  in  Fig,  4.7.  The  gas  leaves  the  cyclone  through 
a  1<y  in,  bore  pipe  having  a  90  lobster-back  bend,  and  exhausts  into  a  six-inch 
bore  pipe  through  a  7°  diffuser. 
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The  lower  end  of  the  14  in.  bore  exhaust  pipe  is  fitted  with  a  22  gauge  skirt 
to  act  as  a  drip-point  within  the  inlet  vortex  for  any  liquid  cesium  film  running 
down  the  outer  wall  of  the  exhaust.  Without  this  skirt,  any  such  cesium  film  would 
drip  off  at  the  exhaust  lip  and  be  carried  away  up  the  tube  by  the  outgoing  helium. 

A  collecting  pot  is  provided  at  the  lower  end  of  the  cyclone  to  receive  the  cesium 
as  it  flows  down  the  core.  This  pet  can  be  drained  into  a  lower  pot  through  a  in. 
bore  pipe  and  a  bellows-sealed  valve. 

All  parts  are  made  of  BN58B  stainless  steel  and  all  Joints  argon  arc  welded. 
Unnecessary  projections  and  roughnesses  inside  the  cyclone  were  removed  to  improve 
the  cyclone  collecting  efficiency. 

if. 2. 3  Pressure  losses 

Before  considering  the  system  pressure  loss  it  is  necessary  to  know  if  the  gas 
Mach  number  is  significant.  The  maximum  average  velocity  will  occur  at  the  cyclcne 
inlet.  The  Mach  number  there  will  be 


Ms 


S 

Aj  V?  P 2  ^2 


Now  G  = 

As  = 
X  - 
P2  = 

P2  = 

so  that  Ms  = 

This  is  sufficiently  low 


6.7  gn/sec 
10.56  cm^ 

1  .63 

0.291  X  1,003  X  10’  X  981 

2.95  X  10’  dyne/cm^ 

3.8  X  10"’  gu/ cm^  , 

O.lifB 

to  be  negligible  in  pressure  loss  calculations. 


The  overall  pressure  change  in  the  separator  system  can  be  conveniently  split 
into  the  following  parts;  inlet  section-changer,  cyclone,  exhaust  pipe  and  bend, 
diffuser.  These  are  considered  separately. 


The  cyclone  inlet  Reynolds  number  is 


where 


Res 
G  = 


As  Ms 

6.7  g^v'^sec 


4.5 


2  B  H 
c  c 


D„  = 


H 


H  +  B 
c  c 


2  X  0.875  X  1>675  x  2,54 
2,75 


Therefore 


=  3 *03  cm 
Aj  =  10,56  cm^ 

^2  =  0,02  5  10“^  poise 

Re 2  =  8,3  X  10^ 


The  inlet  flow  is  thus  turbulent.  The  pressure  loss  in  the  section  changer 
can  therefore  be  expressed®  as 

Pi  -  P2  =TPi2  (V2"  -  Vi"  +  KVj") 


1  G' 


2  Pi  2 


1 

"  Ai 


+  0.05 


l_Gj_  r  1 ,05  _  ”] 

2  Pi  2  _^7  ^  *  J 


(4.3) 


G  =  6,7  gav/cec 
P12  =  3.8  X  10"®  gn/cffi’ 

A2  =  10-56  cm" 

Ai  =  182  cm" 

so  that  the  inlet  section  pressure  loss.  Pi  -  Pj  =  5«58  x  10®  dyne/ cm" 

and  the  cyclone  pressure  loss,  P2  -  P3  =  3A5  x  10^  dyne/ cm" 

_  Jt±. 


8 

T 


and 


In  the  exhaust  pipe  arid  bend  Reynolds  nuxober  R03  =  —  , 

G  =  6,7  fin/ sec 


D  =  4,12  cm 
e 


so  that 


=  0,023  X  10”"  poise, 

4  X  6,7 


Re  3  ~ 


®  "  v  X  L.12  X  0,023  X  10 

=  9.01  X  10® 


FT 


The  pressure  loss  in  the  pipe  is  given  by 


P3-P4  = 


4  f  L 


2  P34  A3 


7  dyne/ cm" 


(4,4) 


V( 

8J 


4,6 


in  which 


f  =  0.008 


5er 


5) 


■) 


=30+32  ^  4.12 
=  161 ,7  cm 
0  =  6.7  gm/sec 


P34  =  3-8  X  10“®  gn/cai^ 

A3  =  13.3  Cffl^  , 

so  that  the  extiaust  pipe  and  head  pressure  loss,  P3-P4  =  4.72  x  IO’  dyne/cm^ 

For  a  diffuser  operating  at  low  Mach  number  the  pressure  change  is  given  by 


P5-P4  = 


(1  -  f )  -  V5" 


dyne/  cm^ 


0^  r  1 

2  P54 


(1  -0 


in  which  P54  =  3.8  x  10“®  gn/cm^ 

G  =  6.7  gn/sec 
A4  =  13.3  CO* 

As  =182  cm* 

C  =  0*088,  (de  Kovats  and  DesouiA  ) 


so  that  the  diffuser  pressure  loss,  P5-'P4  =  3»04  x  IO’  dyne/ cm*  . 


Thus  Pi 
P2 
P3 

P4 

so  that  Pt 


P2  =  +5.58  X  10’ 

P3  =  +3.45  X  10^ 

P4  =  +4.72  X  10’ 

P5  =  -3o04  X  10’ 

P5  =  +41  .76  X  1 0’  dyne/ cm* 

=  0.605  psi 


(4.5) 


This  can  be  expressed  generally  as 

Pi  -  P5  =  1.05  X  10-’  psi  . (4.6) 

4.2.4  Collecting  efficiency 

The  collecting  efficiency  of  a  cyclone  depends  upon  the  cyclone  geometry, 
velocity  of  the  carrier  fluid  at  inlet,  the  nature  of  the  material  to  be  extracted, 
and  the  particle  size.  The  data  used  in  Fig.  4.8  were  abstracted  from  Perry^  and 
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ahowB  the  variation  of  collecting  efficiency  with  a  non->dlm8nsiona?  particle  diameter 
Dp 


for  8  cyclone  as  defined  in  Fig.  4.6.  Dpc  refers  to  xhe  particle  size  for  which 

Considerations  of  the  forces  acting  on  the 


cyclone  has  an  efficiency  of  0,5. 
particles  in  a  cyclone  show  that 


D. 


PC 


w  N  e 


S 

TT 


D 


fPp  “  fil 


(4.7) 


where  Sep  and  Cj^  are  related  by 


8  E'  (pp-p) 

*2''  o  p  n.p* 


•  (4.8) 


For  particle  systems  tiaving  0.0001  <  Re  <  2.0,  Stokes’  law  is  valid,  and  the 
above  equations  reduce  to 

I  rrB  "■ 

®?C  “  2  ®  Ne  V2  (pp  -  p)  . 

Using  the  given  design  conditions,  with  Pig.  4.8,  it  is  possible  to  predict  (see 
table  below)  the  collecting  efficiency  of  the  cyclone  for  cesium  droplets  of  different 
sizes. 


Thus,  since 


then 


fj  =  0,023  X  10"*  poise 

B  !t  2.22  cm 
c 

Ne  =  5  (see  Perry^ ) 

Vj  =  1,67  X  10^  oj^/sec 

Pp  =  1.8  ga/m^ , 

Dp^  =  0.695  X  10"^  OB 


1 

V 

0.5 

0.72 

0.37 

1 .0 

1,44 

0.65 

1.5 

2.16 

0.82 

2.0 

2.88 

0.90 

2.5 

3.6 

0.94 

3.0 

4.32 

0,96 

3.5 

5.03 

0.97 

4.0 

5.75 

0.97 

4.5 

6.47 

0.98 

5.0 

7.20 

0.98 

4.8 
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It  must  be  noted  here  that  these  predicted  collecting  efficiencies  refer  to 
iolid  spherical  particles r  When  the  particles  are  liquid,  the  secondary  vortices 
in  the  upper  annular  section  of  the  cyclone  can  induce  i  flow  of  liquid  along  the 
inner  walls  towards  the  exhaust,  thus  reducing  the  efficiency.  This  effect  can 
be  reduced  by  introducing  drip  points.  If  the  liquid  particles  agglomerate,  then 
the  collecting  e fficiency  will  increase. 

.V.2.5  Pressure  loss  tests 

Before  installation  into  the  MPD  closed-loop,  the  cyclone  separator  was 
tested  in  the  open  circuit  (Fig.  4*9)  for  pressure  loss  and  collecting  efficiency. 

Air  was  drawn  from  the  laboratory  by  a  centrifugal  fan  and  led  by  3  in» 
bcre  piping  into  the  cyclone  inlet.  The  3  in»  pipe  was  extended  into  the  conical 
section  changer  to  ensure  e.  smooth  transition  from  3  in.  diameter.  The  air 
exhausted  into  the  laboratory  through  the  diffuser.  Plow  variation  was  possible  by 
variation  of  the  fan  speed,  or  by  a  gate  valve  located  just  downstream  to  the  fan. 

For  this  test  the  x  in.  bore  pipe  at  the  lower  end  of  the  cyclone  collecting 
pot  was  plugged  with  a  rubber  bung.  The  air  mass  flow  was  measured  with  a  plain 
intake  orifice^  of  outside  diameter  5«72  cm,  attached  to  the  fan  inlet.  A 
paraffin  filled  inclined  manometer  measured  the  pressure  difference  across  the 
orifice. 


The  pressure  difference  across  the  cyclone  system  was  determined  by  measuring 
the  static  pressure  of  the  flowing  air  just  up  tream  to  the  conical  section  changer, 
using  a  water  filled  manometer,  and  assuming  the  exhaust  pressure  to  be  atmospheric 
at  the  diffuser  outlet. 


Measurements  were  taken  for  air  flows  ranging  from  10  to  30  gn/seo,  the 
observed  orifice  heads  and  static  pressures  being  recorded  graphically  in  Fig.  4*10. 
From  these  results 


H  =  36.2  h 


For  the  orifice. 


where  the  units  of  h  are  in.  of  H2O, 

H  are  cm  of  H20« . 

gn/sec  where  the  units  of  h  ere  in.  of  HaO, 

P  are  mm  of  Hg, 

T  are  K. 

ft 


The  system  pressure  loss  is  therefore 


Pi-Ps  =  9.24  X  10" 2  J  dyne/cm* 
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During  the  test 


T  =  296®K 
a 

=  747.65  mffi  Hg, 

so  that  P,  -  P3  =  35.9  G*  dyne/ cm*  (for  air) 

To  determine  the  losses  introduced  by  the  cyclone  itself,  it  is  necessary  to 
separate  the  losses  in  the  section  changer,  exhaust  pipe,  artd  diffuser. 

Prom  equation  (4.3) 


and,  for  air 


so  that 


Pi  2  =  1.17  10"^  ga/om* 

Ai  =  45*7  om^ 

Aj  =  10,56  cm^ , 

Pi-p2  =  3.8  G*  dyne/ cm’  (for  air) 


From  equation  (4.4) 


p  -P  -  tii 

-  D,  2  P34  Aj"  • 


Now 


R6i  s 


Jt± 


V  D,  P3 


therefore 
and  then 


Also 


so  that 


D  >=  4*12  cm 

=  0*0182  X  10"’  poise 
C  =  10  to  30  gn/sec 
Res  =  1 .7  X  10^  to  5.1  X  10^ 
f  =  0.0065  to  0.0053 

Taking  an  average  value 

f  =  0.0059 

P34  =  1.17  X  10"’  gn/cm’ 

As  =  13.3  cm’ 

L  =  161 .7  cm, 

PS-P4  =  2.24  G*  dyne/om’ 


Prom  equation  (4.5) 

2  Ps 


4.10 


so  that 


P54  =  1*17  ><10^  gn^/cffl^ 

A4  =  13.3  cm^ 

A5  =  182  cm^ 

C  =  0.088, 

P5 -P4  =  2,57  dyne/ cm* 

Now  the  cyclone  pressure  loss  is  (Pa-Ps)  and 

(Ps-Pj)  =  (P1-P3)  “  (P1-P2)  -  (P3-P4)  “  (?4-P5) 

=  (35.9  -  3.8  -  2.24  +  2.57)  dynes/cm^ 

(P2-P3)  =  32.43  G*  dynes,/ cm* . 

Expressed  in  inlet  velocity  heads, 

(P2-P3)  =  8.45  hj. 

4.2.6  Collection  efficiency  testa 

To  simplify  the  experimental  teciiniques,  che  cyclone  was  tested  with  a  mixtiire 
of  air  and  carhonyl  iron  powder  simula  ting  cesium  (manufactured  by  Inco  Mond) .  the 
particles  being  spherical  and  of  known  sizes.  Two  grades  of  powder  were  used; 

MDG  and  MCHP  of  average  particle  sizes  2.9  and  7  micron. 

The  iron  powder  was  injected  into  the  air  through  the  static  pressure  tapping 
upstream  to  the  cyclone  section  changer  in  the  following  maiuier.  Powder  was 
placed  in  a  pot  fitted  with  a  perspex  lid  and,  by  blowing  dried  nitrogen  through  a 
1/16  in.  bore  pipe  tangentially  into  the  pot  base,  the  powder  could  be  fluidized  and 
car^*iod  through  an  4  in*  boro  pipe  into  the  main  air  flow.  To  collect  the  powder 
separated  by  the  cyclone,  a  weighed  glass  boiling  tube  was  attached  to  the  4  in. 
pipe  at  the  low^r  end  of  the  cyclone  by  a  bored  rubber  bung. 

During  a  test  run,  a  known  weight  of  either  MD&  or  MCHP  grade  powder  was 
placed  in  tho  clean  fluidizor  pot  and  the  lid  sealed.  The  main  air  flow  was 
adjusted  to  tho  required  value.  Then  tho  nitrogen  flow  into  the  fluidizsr  was 
started  and  adjusted  until  tho  powder  was  seen  to  be  thoroughly  entrained  in  the 
induced  vortex.  Shortly  after  this  condition  was  achieved,  powder  could  bo  saen 
Settling  at  the  bottan  of  the  glass  collecting  tuba.  Prom  time  to  time  during  tho 
run,  tho  cyclone  was  struck  with  a  plastic  hammer  to  release  any  powder  clinging 
to  the  cyclone  inner  walls.  When  the  r*  xulzing  pot  was  empty  the  air  flow  was 
increased  to  a  maximuTn  to  clean  the  cyclone  of  collected  powder.  Tho  collecting 
tube  was  detached,  closed  with  a  weighed  bung,  and  weighed.  This  procedure  was 
repeated  for  various  air  flows  using  each  of  tho  powders. 
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After  completion  of  the  moesurementa  the  cyclone  and  section  changer  were 
carefully  removed  from  the  circuit.  Iron  powder  had  settled  on  the  bottom  of  the 
pipe  belo^  and  downstream  of  the  injection  point  in  the  upper  wall  of  the  pipe. 
Powder  had  also  collected  in  the  dead  space  between  the  section  changer  none  and 
the  length  of  >  in.  pipe  protruding  into  the  section  changer  (Pig.  4.9)*  The 
powder  was  deposited  in  layers  of  light  grey  and  brown  grey,  showing  that  both 
grades  were  present.  All  of  this  powder  was  collected  and  weighed. 

The  measured  air  flows,  injected  weights  of  powder,  and  collected  weights 
of  powder  are  recorded  in  the  following  table. 


gn/sec 

Powder 

grade 

®P  > 
micron 

mm 

n 

®PC’ 

micron 

■ETn 

^PC 

1 

9.1 

MCMD 

V 

110 

31.2 

78.8 

78.8 

0.997 

1.52 

4*6 

2 

13*9 

UDfr 

2.9 

23.1 

0.008 

93.1 

22.4 

0.97 

1.25 

2.32 

3 

32.5 

MCHD 

7 

69.2 

5.38 

63.8 

75.4 

1.15 

0.806 

00 

00 

9.2 

MDG 

2.9 

79.5 

4.2 

75.3 

51.5 

0.686 

1 .46 

1.98 

5 

5.65 

MCHD 

2.9 

19.5 

1.69 

17.8 

14.1 

0,794 

1.935 

1.5 

6 

9.32 

MCHD 

7 

48.2 

13.7 

34.5 

31.2 

0.92 

1.415 

7  95 

7 

33.3 

MDG 

2.9 

23.2 

0.34 

22.9 

23.4 

1.02 

0.7.7 

' »'. 

8 

5.65 

MCHD 

7 

36.^ 

16.7 

19.2 

21 . 

1 .09 

1.935 

3.62 

Before  calculating  collecting  efficiencies  it  was  necessarj*  to  correct  the  measured 
injected  weights  for  the  amount  of  powder  lost  by  settling  dvorim  ach  run.  This 
was  done  by  proportioning  the  total  deposited  weight  according  to  the  efficiency 
as  a  gravity  settler  of  the  3  in*  pipe  downstream  to  the  injection  point  for  each 
set  of  test  conditions.  According  to  Perry*,  a  rectangular  gravity  settler  has 
the  efficiency;  ^ 

t 

""  n  Y  * 
s 


The  weight  of  powder  lost  by  settling  during  each  test  was  thus  estimated  as 


W 

8 


*1 

2  ' 


2  W 

8 


The  weight  of  powder  actually  entering  the  cyclone  was  calculated  as 


and  the  collecting  efficiency  as 

W 


c 


4.12 


The  ratio  ^P/D^.  was  evaluated  for  each  test  using  equation  (4.9).  The 
PO 

appropriate  values  of  and  Rsp  were  found  from  equation  (4.8),  in  c'^njunction 
with  Fig.  114  of  reference  4  which  shows  the  variations  of  drag  coefficient  with 
Reynolds  number  for  spheres. 

The  results  of  these  tests  are  presented  in  the  above  table  and  in  Fig.  4.11. 
The  ,  ^P/Dp^  characteristic  recommended  by  Perry^  is  also  given.  The  two  curves 
do  not  differ  significantly,  so  that  the  design  is  s'^tlsfactoryw 

On  the  basis  of  the  measured  characteristic,  the  expected  collecting 
efficiency  has  been  calculated  for  the  design  conditions  <^Fig.  4.12). 

4.2.7  Loop  operation 

The  unit  was  fitted  into  the  loop  (Fig,  4.5)  at  the  heat  exchanger  outlet;  the 
perfcrmi.anoe  from  seeded  runs  RHT  35  and  55  (and  later  runs  not  reported)  showed  ein 
actual  esaium  removal  efficiency  of  -'•12  per  cent.  A  commercial  Centrifix  Separator 
has  been  ordered  and  will  be  installed  in  the  loop  for  seeded  runs  planned  early  in 
December  1964.  A  removal  efficiency  better  than  99  per  cent  is  guaranteed. 

4.5  ELECTROSTATIC  PRSCIPrrATION 

Careful  consideration  has  been  given  to  the  feasibility  of  a  cesium  remcval 
system  employing  electrostatic  precipitation  techniques  for  the  oesium-helium  mist. 
Such  a  system  with  negligible  pressure  loss  and  very  high  aved.lable  efficiency  is 
extremely  attractive  for  application  in  the  closed  cycle  KPD  power  generation 
facility.  Whereas  with  a  cyclone  separator  the  separation  efficiency  for  gas-bome 
particles  of  less  than  SO  micron  falls  off  rapidly,  the  electrostatic  precipitator 
maintains  a  high  separation  efficiency  for  small  particles. 

A  theoretical  feasibility  stud^-^  has  been  carried  out  and  small  development 
units  built  and  tested.  At  present , stable  coronas  have  been  achieved  with  helium 
alone  (Fig.  4.5)  but  insufficient  operation  has  been  possible  with  cesium-helium 
flows  to  prove  performance  potential, 

4.5*1  Principles  and  theory 

In  electrostatic  precipitation  unwanted  charged  particles  drift  out  of  the 
main  gas  flow  under  the  influence  of  an  electric  field,  the  efficiency  of  removal 
depending  upon  the  applied  voltage.  Particles  must  be  electrically  charged,  one 
means  of  achieving  this  being  to  create  a  corona  by  applying  a  sufficiently  strong 
field  between  the  electrodes.  Two  types  of  d.c,  corona  can  be  set  up,  depending 
upon  whether  the  corona  electrode  is  positive  or  negative.  Positive  corona  is 
sot  up  at  higher  voltages  than  for  negative  corona  in  hydrogen,  carbon  dioxide  and 
nitrogen;  in  helium,  oxygen  and  air  the  critical  eleotrio  field  is  lower  for  positive 
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than  for  negative  corona  at  pressures  greater  than  1  era  Hg,  but  higher  at  lower 
pressures.  However,  positive  corona  is  usually  employed  because  of  its  mere 
uu^fona  distribution;  negative  corona  appears  as  tufts  or  beads  along  the  corona 
electrode,  the  presence  of  which  complicates  txie  precipitation  mechanism  and  renders 
theorst:?cal  analysis  difficult. 


Charging  of  the  particle  takes  place  by  ionic  impact,  either  by  attachment  to 
gas  ions  migrating  to  the  cathode  or  with  ions  under  thermal  motion  (effective  only 
for  particles  of  less  than  x  micron  diameter). 


A  full  analysis  shows  that  the  number  (n)  of  elementary  charges  (e)  acquired 
in  time  t  by  a  spherical  particle  under  migrating  ion  impact  conditions  is  given  by 


n  = 


n  E  a' 
"  c 


t  + 


(4.10) 


wNek 


whore  p  is  a  constant  depending  upon  the  particle  material, 
a  is  the  radius  of  the  particle, 

is  the  field  strength  at  corona  olo  trode  surface, 
N  is  the  ion  concentration,  surd 
k  is  the  ion  mobility. 

As  t  -*  09,  then  - 1 


In  practice  maxtmam  charge  is  attained  in  a  fraction  of  a  second,  and 

pB  a* 

“(MAX)  =  *■  V“ 

or  maximum  charge  q  = 


so  that  q  =  pEa*.  . (4.1 1) 

c 

The  constant  p  =  3  for  conducting  pairticlos  and,  for  non-conducting  particles 
p  =  where  c  is  the  dielectric  constant. 


Under  conditions  of  thermal  ion  motion,  the  number  of  charges  accepted  by  a 
particle  is 


a  KT  ,  /  ^  IT  as  Ne‘t\ 


(4.12) 


where  K  is  Boltzmans  constant, 

T  is  absolute  temp, 
s  is  r.m.s  ion  velocity,  and 
N  is  ion  concentration. 
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For  concentric  cyliniers  it  has  been  shoTO  that  a  corona  can  only  be 
produced  when  the  ratio  of  outer  to  inner  radii  of  electrodes  exceeds  a  certain 
critical  limit;  the  starting  voltage  for  a  corona  depends  upon  the  si&e  and  shape 
of  the  corona  electrode  (wire)  anri  interspace. 


The  field  strength  B  between  two  radii  and  can  be  defined  by 

r 


V  = 


B  dr 


R; 


where  V  is  the  potential  difference  between  E^  and  Rj . 
The  field  strength  at  a  radius  r  is  given  by 

V 


E  = 


r  log  R2/R1 

O 


c 

r 


(4.13) 


where  c  =  7*-= - J  constant 

logg  R2/R1 

In  the  presence  of  an  ionic  apace  charge  the  field  is  modified  to 

^  +  K  -  k-va  =  0  .....  (4.14) 

dr  r 

where  a  is  space  charge  per  unit  volume. 

At  radius  r  the  ionic  current  i  per  unit  length  of  electrode  wire  is  given 
by 

i  =  2  w  r  <r.  (ion  velocity) 

=  2  F  r  (T  EX  . (4.15) 

where  k  is  ion  mobility. 

Prom  equations  (4.14)  and  (i;.15)  the  following  equation  may  bo  obtained 

,  /2i  o2^4 

B  =  (y  -  it)  ’ 

and,  for  largo  values  of  i  and  r,  —y  can  be  neglected 

so  that  S  =  (y)"^  .  (4.16) 

When  particles  are  present  equation  (4.16)  is  modified  to 

B  =  (^)  ^  .  (1  +  J  P  S  r)  . (4.17) 


where  S  is  the  surface  area  of  particles  per  unit  volume. 

This  equation  shows  a  linear  relationship  between  field  strength  and  radius,  which 
is  dependent  upon  the  nature  and  concentration  of  particles.  Pig.  4.13  shows  how 
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particles  effect  the  field  strength. 


The  force  acting  on  a  particle  with  charge  q  in  a  field  is  given  by 


Fran  equation 


F.  =  p  K  S  a* 
c  P 


where  a  for  a  high  particle  concentration. 

The  viscous  drag,  by  Stoke" '  law,  is 

Fa  =  6  va  p  vT  w 


where  a  is  particZ  i  radius , 

p  is  viscosity  of  gas,  and 
w  is  velocity. 

The  transverse  terminal  velocity  u  is  reached  fairly  quickly,  so  that 


Fi 

Therefore  u 

For  a  conducting  particle 

u 

The  longitudinal,  vslocity 

V 


=  P2 

p  S  £  a 
6  tr  p 


p  =  3  and  therefore 
a  S  E 

=  0.16  — ^ 


for  ■’aminar  flow  v  is  given  by 
=  2V  (1  -  ) 


(4.18) 


(4.19) 


where  V  =  average  flow  of  gas 

y  =  radius  from  centre 
R  =  outer  radius  of  tube, 

A  particle  clone  to  the  inner  electrode  will  take  the  longest  time  cf  all 
to  precipitate,  h-u"'  will  give  a  theoretical  maximum  for  the  lerigth  of  ttiho  required. 

The  longitudinal  velocity  of  such  a  particle  =  2V 

at  y  =  0  . (4-20) 

and  xhe  length  of  tuba  necessary  to  precipitate  a  particle  is  given  by 

^  -  3  •  u 

and  from  eqxiations(4.l8),  (4.20)the  maximum  length  of  tube  required  la  given  by 

r  _  it  2  V  R  p 

(MAJC)  -  3  •  0.16  r 
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where 


is  tha  minimum  radius  of  particle  in  the  system. 


Therefore 


T  16.7  V  R  u 


(4.21 ) 


If  V  is  in  cn/sec,  a  and  in  cm,  ^  in  poise, 
(300  volt/c®  s  1  esu/cm)  then  L  will  be  in  cm. 


E  and  2  in  osiV^cm 
c  P 


Sample  Calculation 

Without  full  investigation,  both  theoretical  and  experimental  >  the 
critical  voltage  at  which  the  corona  appeal’s  cannot  be  determined.  The 
potential  difference  to  precipitate  all  cesium  particles  can  bo  calculated  from 
the  following  equation. 

T  16.7  V  : 

^(MAX)  “  Eq  *p 

The  following  dimanai  nr.s  and  flow  data  are  used: 


=  2  m 

=  7  cm 
=  4.. 5  nv^aoc 
=  25*5  X  10”®  poise 
=  1  Micron 

The  product  3  is  found  to  bo  67O  (esH/cm)^,  If  it  is  assumed  that  2  «  K 
1  P  c  p 

then  E  «(67C)‘  »  26  osVcm  =  7800  volt/ cm  and  the  potential  difference  required 
c 

is  about  40  kxlovolts. 

The  above  calculation  allows  for  extreme  conditions;  an  increase  in  minimum 
particle  size,  or  a  decrease  in  temperature  "  nd  hence  viscosily),  flow  rate,  or 
radius  of  the  tube  would  lower  the  potential  difference  needed  to  produce  the  corona. 

4.4  ULTRANSONIC  CSSITJM  REMOVAL  SYSTEMS 

Ultrasonic  techniques  are  also  being  considered  to  promote  cesium  di’oplet 
agglomeration  and  facilitate  its  removal  from  the  helium  flow. 


maximum  length  for  precipitation 


reuiius  of  tube  R 

mean  helium  flow  V 

viscosity  of  helium  at  150°C  ^ 

radius  of  smallest  farticlo  ®(MIN) 


For  the  next  series  of  seeded  runs  a  500  watt  power  source  is  to  bo  used  to 
drive  an  ultrasonic  transducer  titanium  31  8A  alloy  probe  (Fig.  4.4) »  Tx*anamitted 
powers  of  up  to  1  watt  are  expected.  The  pmbe  is  to  be  mounted  at  the  bottom 
of  a  vertical  glass  column  (6  in.  diameter,  8  ft  high)  and  is  arranged  to 
propagate  axislly  in  the  direction  of  the  helium-cesium  flow.  Stringent 
pre-seeding  tests  will  be  carried  out  to  determine  possible  fatigue  characteristics 
of  the  glass  system. 
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Should  tho  electrostatic  precipitator  and  ultrasonic  systems  prove  practically 
feasible,  a  combination  of  the  two  will  give  extremely  high  efficiency  of  seed 
removal . 

4.5  CSSIUM  VAPORIZATION  AND  SUPERiSAT 

Until  the  problem  of  seed  removal  is  resolved,  only  short  duration  seeded 
injection  will  be  possible.  For  the  immediate  future  a  considerable  programme  is 
planned  for  investigations  at  nozzle  inlet  temperatures  lower  than  1800^K.  A 
system  has  been  designed  (Fig.  4.5)  (and  is  in  the  final  development  stage)  to 
temporarily  replace  the  high  temperature  helivua  heater,  high  temperature  transfer 
«ind  cesium  vaporizer.  The  system  is  horizontal,  axially  in  line  witn  and  connected 
directly  to  a  subsonic  nozzle  and  generator  duct.  It  comprises  a  stainless  steel 
mesh  electrical  heater  coruiected  mechanically  in  series  with  a  boron  nitride 
insulated  tungsten  element  electricaJ  resistance  heater,  the  tungsten  element  being 
in  three  stages  to  allow  for  power  balance.  Cesium  is  to  be  injected  through  three 
Jets  and  atomized  with  the  helium  flow  prior  to  entry  into  the  stainless  steel  mesh 
heater  (available  power  27  kilowatts) which  takes  helium-cesium  to  “900°C  (cesium 
in  vapour  state).  The  second-stage  tungsten  heater  (power  20  kilowatts)  raises 
the  helium-cesium  temperature  (mass  flow  helium  5  gni/sec  at  1  ata;  cesium  fraction 

1.5  ^/o)  to  approximately  1800°K. 

Careful  consideration  has  been  given  in  the  design  to  the  problems  of  electrical 
breakdown  due  to  the  presence  of  cesium.  Tho  unit  is  to  be  commissioned  shortly 
and  if  satisfactory  should  give  additional  data  cn  the  problems  of  duct  surface 
contaminat-’  on. 

R]£FEm?C1SS 

1  P4NKHURST,  R.C.  and  HOLDER,  D.W.  Wind  tunnel  technique.  Sir  Isaac  Pitman 
and  Sons  Ltd.  1952 

2  McNAB,  IJ".  and  COOPER,  N.A.  Flow  processes  in  MPD  generators.  IRD  65-82 

1965 

5  He-He  recuperative  heat  exchanger  for  the  IRD  MPD  circuit.  IRD  Tech,  Memo, 
ref.  TA/MJB,  20th  June  1965 

4  PERRY,  J.H.  (Editor}.  Chemical  Engineers*  Handbook.  5x*d  ed.  1950.  McGraw- 
Hill  Book  Co.  Inc. 

5  McADAMS,  W.H.  Heat  transmission.  ed.  1954.  McGraw-Hill  Book  Co.  In'', 

6  DE  KOVATS,  A.  and  DESMUR,  G.  Pumps,  fans,  and  compressors.  Blaokie  1958 

7  LINFORD,  A.  Flow  measurement  and  meters,  2nd  ed.  E.  and  F.I',  Spon  Ltd, 

1961 . 


4.18 


NOm^CLATDRE  (Chapter  k) 


A  fluid  flow  area 

B  magnetic  field 

cyclone  inlet  width 
Cjj  drag  coefficient 

cyclone  inner  diameter 
cyclone  exhaust  diameter 
Dp  particle  diameter 

Dp^  cut  particle  diameter 

F  frictional  losses 

f  fanning  friction  factor 

G  mass  flow 

g  gravitational  acceleration 

h2  cyclone  inlet  velocity  head 

h  orifice  pressure  difference 

H  manometer  head 

cyclone  inlet  height 

height  of  gravity  settler 

K  dimensionless  coefficient 

L  pipe  length  for  friction  loss 

length  of  straight  pipe 

Lp  equivalent  length  of  bend 

L  length  of  gravity  settler 

s 

K  Mach  number 

Ne  No.  of  turns  by  particle  in  cyclone 

P  pressure 

Re  Reynolds  number 

r^  radius  of  particle  path 

T  temperature 

u^  terminal  velocity  of  particle 

V  velocity  of  gas 

injected  weight  of  powder 
W_,'  corrected  injected  weight 

settled  weight  of  powder 
collected  weight  of  powder 
^  ratio  of  specific  heats 

L  differential 

t  diffuser  efficiency 

rj  collecting  efficiency 

77  gravity  settler  efficiency 

8 


/i  viscosity 

f  diffuser  losses 

p  density 

Pp  particle  density 

Suffices 

0  refers  to  generator  duct  outlet 

1  refers  to  heat  exchanger  inlet 

2  refers  to  cyclone  rectarjtular 

inlet 

3  refers  to  cyclone  outlet 

4  refers  to  exhaust  diffuser  inlet 

5  refers  to  exhaust  diffuser  outlet 

P  refers  to  particle 

f  refers  to  flui^ 

a  refers  to  ambi  :,t  conditions 
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CESIUM  VAPORIZER 


ELECTROSTATIC  PRECIPITATOR  AND  CYCLONE  SEPARATOR 


FIG  4-3 


^il«  Hi  ' 


A 


ULTRASONIC  TRANSDUCER  WITH  TITANIUM  PROBE 


Arbitrary 


CYCLONE  TEST  CRCUfT 


FIG  4 .9 


COCUECTOR  TUBE 


Bouir'i  y 


CYCLONE  PRESSURE  LOSS  TEST  RESULTS 


FIG  410 


CYCUOKtC  P»B?Ca»U«fc  U05ts,  CM 


■■•ItllilllUtSI,' 


FIG  41 3 


EFFECT  OF  OUST  BURDEN  ON  FIELD  STRENGTH 


CHAPTER  5 


GENERATOR  DUCT 


by 

R*Browii,  W^JpSlator,  V.Brcswn,  K.Casa,  J, Davidson  an^  G. Johnson 
5*1  SPECIFICATION 

The  generator  duct  system  employed  in  the  cesium-seeded  runs  in  March  196^ 

(RHT  52,  see  TRD  64.-36,  Bibliography  number  45)  had  exhibited  significant  electrical 
end  leakage  to  earth.  For  the  hot  run  (RifT  5^)  and  the  subsequent  cesium-seeded 
runs  (RHT  35)  (see  Chapter  3)  the  generator  duct  walls  were  extended  axially  at 
the  subsonic  diffuser  end  and  an  electrically-insulating  expansion  nozzle  in  boron 
nitride  was  fitted.  Fig,  5»1  shows  the  tantalum  (0.02  in.  thick)  high  t^nperature 
transfer  section  and  cesium  vapourizer,  with  the  attached  boron  nitride  nozzle. 

Duct  walls 
Duct  dimensions 

Electrodes 

Probes 

Expansion  nozzle 
Thermal  insulation 
Thermocouples 

Fig.  5.3  shows  the  assembly  (diffuser  end)  loaded  into  lha  generator  case. 

Fig,  5.4a  shows  the  generator  assembly  after  RHT  34  ard  35  in  a  partly  dissraabled 
state.  Bonding  of  the  boron  nitride  plates  had  occurred  at  the  clo.se- toleranced 
interface  joints;  this  was  attributed  to  ingress  and  reaction  of  Refrasil  dust. 

The  plates  were  finally  forced  apart  mechanically. 


boron  nitride 

6  in.  X  l4  in.  x  ^  in.  (^-t  in.  diverging 
diffuser) 

16  pairs  ^  in.  semi-circular  section  tungsten 
10^  tantalum  90^  alloy  (Fig.  5.2a) 

~  in.  diameter  tantalum,  two  on  each  side  of 
nozzle  throat  and  four  at  diffuser  end  of 
generator 

boron  n.tride  profiled  for  Mach  2.3  operation 
(Fig.  5.2b) 

generator  duct  and  transfer  section  Zirconia  ’E‘ 
fibre  with  outer  wrappings  of  'Refrasil’  tape 
in.  tantalum-sheathed,  tungsten-tungsten  26^ 
rhenium  '  upstream  of  nozzle  in  gas  stream,  2 
flush  with  duct  wall  downstream  of  generator 
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5.2  KLECTRICALLY -CONDUCTING  LAYKRS 

A  fairly  uniform  electrically-conducting  layer  (less  than  1  ohcv/cm)  was 
observed  on  all  four  duet  inner  wall  surfaces  and  the  inferior  surfaces  of  the 
expansion  nozzle  (Figs.  5.5a  and  5.5^).  This  phenomenon  has  been  discussed 
previously  (iRD  64-36)  emd  extensive  analysis  of  electrically  conducting  surface 
layers  has  been  carried  out  without  determining  the  mechanism  of  layer  transport 
or  the  major  components  contributory  to  the  electrically  conducting  properties  of 
the  layer. 

Tests  have  been  conducted  on  small  samples  of  electrically-conducting 
generator  walls  in  both  boron  nitride  and  recrystallized  alumina,  heating  for  one 
hour  under  controlled  atmospheres  from  vacuum  (<10""^  torr)  to  1  ata  helium  pressure, 
through  a  temperature  range  1450°K  to  1850°K.  Samples  ta^en  from  ducts  used  in 
runs  before  RHT  33  have  lost  the  conducting  layer  with  this  procedure  and  recovered 
the  initial  electrical  resistivity.  A  further  test  included  fitting  six  electrically 
conducting  samples  of  recrystallized  alumina,  three  upstream  and  three  downstream 
of  the  cesixin  injection  station  in  the  MPD  tantalum  high  temperature  transfer  section 
allowing  the  layered  surface  to  be  'washed*  by  the  hot  helium  and  helium-cesiim 
flow  during  cesium-seeded  RHl'  36  (not  reported  here).  Inspection  of  the  samples 
after  the  rxxns  showed  that  the  conducting  layer  had  been  removed  in  all  cases  with 
recovery  of  the  surface  electrical  resistivity  (Fig.  5*7).  The  generator  walls 
(boron  nitride)  from  RHT  33»  34  and  35  had  become  layered  and  electrically 
conductive  (less  than  1  ohn/cm).  Prom  these  tests  it  was  deduced  that  the 
conducting  layer  was  unstable  above  1450°K;  it  therefore  appears  that  the  layer 
formation  is  by  deposition  on  the  cooler  generator  duct  walls. 

5.3  ELECTRODES 

Inspection  of  the  tantalum- tungsten  strip  electrodes  after  RHT  34  and  35 
showed  extensive  transverse  creeping  (Fig.  5*6),  Complete  fractures  were  noted 
in  11  of  the  32  strip  electrodes  ^nd  appeared  to  occur  in  each  case  at  the  lead-out 
bend  where  the  internal  stress  promoted  through  forming  would  be  highest.  Some 
inconsistencies  of  power  generation  results  (Chapter  3)  along  the  axis  of  the 
generator  can  possibly  be  attributed  to  the  high  interface  resistance  of  the 
fractures  and  in  some  cases  intermittent  open  circuit. 

A  transverse  section  was  taken  from  a  cracked  electrode  at  a  position  midway 
along  the  acti  •©  electrode  and  the  areas  were  examined  metallographically. .  A^ 
the  electrode-gas  interface  the  layer  had  an  even  thickness  of  17  microns.  At 
the  boron  nitride-electxx)de  interface,  although  there  was  no  evidence  of  a  definite 
layer,  the  electrode  surface  was  extensively  cracked  both  axially  and  transversely 
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(attributed  to  reaction  between  the  electrode  and  boron  nitride).  At  the  comer 
connecting  the  boron  nitride-electrode  interface  to  the  gas-electrode  interface, 
the  depth  of  penetration  of  cracks  was  about  250  microns.  All  three  samples 
wore  extremely  brittlt-. 

An  overall  X-ray  analysis  of  the  electrode  surfaces  showed  the  presence  of 
TajOs  and  TaB.  The  layer  at  the  boron  nitride-electrode  interface  was  assumed 
to  be  TaB,  Samples  of  boron  nitride  and  boron  nitride  surface  layers  were 
malysed  spec  'rographioally.  Unused  samples  of  boron  nitride  showed  -mjcr  R, 
trace  Si  and  minor  traces  A1  a'd  possibly  Cu.  Boron  nitride  from  the  duot  wall 
'*emote  from  the  gas  interface  (non-conducting  electrically)  showed  the  same 
analys:'  i , 

At  the  interface  between  boron  nitride  and  the  eleotrioally"  conducting  layer, 
the  analysis  gave  major  B,  trace  A1  (higher  than  previous  samples),  Ca,  auid  Si  and 
minor  traoes  Pe,  Cr  and  possibly  Cu. 

5.4  &B?IERAT0R  DUCT  WALL  HEATBR  PgVBLOPICNr 

Although  efficient  insulating  materials  have  been  used  externally  on  the 
generator  section  (for  example,  Ziroonia  'S'  fibre)  there  is  a  space  limitation 
transversely  between  the  poles  of  the  magnet  and  it  has  been  recognised  that  the 
generator  walls  (and  especially  the  electrodes)  have  been  at  significantly  lower 
temperatures  than  the  plasma  in  the  generator.  This  has  resulted  in: 

eleotrlcally-conduoting  layers  on  the  generator  wall  inner 
surfaces,  causing  low  inter-electrode  resistance  and  effective 
short  oirouiting  of  the  generator  with  progressive  and  soaetimes 
rapid  performance  deterioration  and  unoartainty  of  results  during 
seeded  runs; 

limitation  of  electron  astission  from  the  electrode  surfaces;  and 

uncertainties  of  gas-wall  interfaoial  temperatures  and  subsequent 
difficulty  in  theontioel  interpretation  of  results. 

For  these  reasons  a  programme  has  been  initiated  for  the  developsent  of 
oxternally»*heated  generator  duot  walls.  Design  considerations  have  been  for  a 
minimum  useful  heater  Dife  of  25  hr  at  up  to  1900°K,  unlfoxmi  temperature 
distribution,  erd  provision  for  lead-through  of  the  required  complex  of  electrode, 
pressure  euid  thermocouple  station  leads.  Barly  work  was  carried  out  using  a  boron 
nitride  duct  top  wall  (7r  in.  x  1-y  in.  x  4  in.)  split  and  screwed  to  accommodate 
single  refractory  metal  (0»010  to  0*020  in.  W  and  Ta)  heating  elaments  (Figs.  5»8a 
and  5*8b). 
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A  rectangular  duct  ha?  been  built  to  the  cross  sectional  dimensions  of  the 
system  used  in  RHT  54  and  35 »  but  with  the  length  reduced  to  4  in.  Internal 
tantalum  (0.020  in.  thick)  heater^  were  fitted  in  each  wall  and  series-connected. 
The  system  was  Insulated  externally  with  zirconia  fibre  and  themally  cycled  to 
1750°E  under  high  vacuum  conditionsj  the  temperature  was  held  a'.  1750°K  for  25  hr 
(Pigs.  b  and  c).  Inspection  showed  no  deterioration  of  the  system  other 

than  embrittlement  of  the  element  due  to  reaction  with  the  boron  nitride. 

The  tests  are  now  being  extended  to  s  tudy  generator  wall  heater  performance 
under  conditions  of  helium  flow  and  at  the  time  of  writing  a  complete  heated 
generator  system  is  being  built  for  power  generation  runs  in  the  closed  loop. 


BORON  NITRIDE  SUPERSONIC  NOZZLE 


1 


FIG  52  b 


GENERATOR  ASSEMBLY  AFTER  RUN 


FIG  5  4a 


HELIUM  LEAKAGE  PATH  BETWEEN  DUCT  END  AND  DIFRJSER 


FIG  5-4  b 


SPLIT  GENERATOR  DUCT  SHOWMG  ELECTRICALLY  CONDUCTING  LAYERS  AND 
PITTING  OF  BORON  NITRIDE  SURFACE 


FIG  5-5a 


A  \ 


SPLIT  BORON  NITRIDE  NOZZLE  SHOWING  ELECTRICALLY  CONDUCTING  LAYERS 

FIG  5-5  b 


CRACKING  OF  90 °4  TANTALUM  KD  4  TUNGSTEN  ELECTRODES 


FIG  56 


ALUMINA  SAMPLES  SHOWING  REMOVAL  OF  ELECTRICALLY  CONDUCTING  LAYER 
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DUCT  WALL  HEATER  DEVELOPMENT  AT  !400°C 

FIG5-8  b 


DUCT  or-  25  hr  L>fE  TEST  AT  I500®C  WALL  TEMPlR/aTURE 

F!G  5-9  b 


DUCT  STRIPPED  SHC.VING  ONE  WALL  HEATER 
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FIG  5-9  c 
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CHAPTER  6 


INSTRUMENTATION 


by 

P.3.  Jones,  M.  Watson,  I.  Wilson,  R.  Brown  and  I.R.  McNab 


6.1  INTRODUCTION 

As  indicated  in  the  previous  Technical  Summary  Report  (IRD  64-56,  Chapter  10) 
much  effort  has  been  devoted  to  accurate  measurement  of  the  plasma  flow  parameters 
in  the  MPD  generator  duct  itself.  Of  particular  importance  is  the  problem  of 
temperatui'e  measurement  for  which  no  satisfact  .'ly  solution  has  been  resolved 
although  the  method  described  in  Ssction  6.4  is  of  some  interest;  consequently  for 
this  and  other  reasons,  heated  duct  walls  to  be  used  so  as  to  operate  under  as 
nearly  isothermal  conditions  as  possible. 

The  system  for  generator  output  parameters  (described  in  Section  6.2)  has 
proved  satisfactory  and  is  currently  being  modified  for  more  sophisticated 
measurements  of  generator  characteristics. 

6.2  GENERATOR  CHARACTERISTICS 

During  cesium-seeded  operation  it  is  essential  to  extract  the  maximum 
information  in  the  available  time.  Cesium  may  only  be  injected  for  5  *«c  periods, 
and  in  the  early  runs,  the  instantaneous  concentration  fluctuated  rapidly;  ideally, 
the  MPD  generator  loading  would  be  cycled  faster  than  these  fluctuations. 

The  following  information  was  required; 

V-I  characteristic  of  each  electrode  pair  under  power  generation 
conditions; 

voltages  on  5  probe  pairs  outside  main  MPD  duct; 

Hall-effect  emf's  and  short  circuit  curTsnts; 
interelectrodte  leakage  resistance; 
plasma  passive  V-I  characteristics; 

electrode  resistance  as  an  index  of  electrode  t''jperature  distribution 
along  the  duct; 

cesium  injection  pert.od;  and 

flow  of  cesium  at  heat  exchi^nger  exit. 
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6*2.1  ¥-1  power  characteriatica 

A  typical  l8-©lectrode  pair  (A,  B,  C....R)  channel  (designated  X)  is  shown 
at  the  top  of  the  circuit  (Fig.  6,1 ).  The  electrodes  (shown  diagrtuamatically  on 
the  extreme  top  left  are  connected  to  the  temperature-scanning  uniselector  and 
current  loop  (described  later)  and  one  lead  from  each  electrode  is  taken  to  the 
patch  panel.  This  has  leads  fitted  with  stacking  plugs  to  allow  parallel 
connection  (any  plug  will  fit  into  the  top  of  any  other);  at  present,  electrode 
pairs  are  each  connected  to  one  loading  and  measuring  channel  by  linking  X11  to 
X13  and  X12  to  X14. 

This  links  the  electrode  pairs  to  a  level  on  one  of  two  25“way  11 -level 
uniselectors  (USA  and  USB)  which,  on  the  first  five  positions,  step  through  five 
load  resistors  (R6  -  R10)  repeated  on  all  18  channels.  At  this  stage  the  groups 

of  relays  (A  and  B)  are  not  energised  and  a  galvanometer  in  the  25“channel  recorder, 
connected  as  a  voltmeter,  is  connected  across  the  electrode  pair.  As  in  Fig,  6,2a, 
a  shunt  resistor  (R4  plus  R5)  is  provided,  also  a  aeries  resistor  selected  from  the 
chain  R1-3  by  a  range  switch.  The  diodes  across  the  galvanometer  have  a  high 
resistance  at  0.1  V,  but  would  shunt  the  galvanometer  and  protect  it  under  overload. 
On  their  sixth  step  the  uniselectors  disconnect  all  loads  so  that  open  circuit 
voltage  is  recorded  in  all  channels. 

Thus  for  each  electrode  pair  the  voltage  at  various  known  loads  is  recorded, 
the  current  in  each  load  can  be  calculated  and  the  V-I  characteristic  derived, 

6.2.2  Probe  voltages 

The  probe  pairs  outside  the  main  duct  are  connected  to  galvanometers  in  the 
6-channel  recorder  with  range  resistors  and  shunts  as  in  Fig,  6.2a.  The  open 
circuit  voltage  is  monitored  continuously  in  this  way;  no  loads  are  connected  to 
the  probes. 

6.2.3  Hall  effect 

On  the  seventh  and  eighth  steps  the  electrodes  of  each  pair  are  shorted  to 
each  other.  Alternatively,  by  moving  flying  leads,  any  of  the  load  resistors  can 
be  oonrected  in  at  this  stage.  On  the  seventh  stop,  the  B  group  of  relays  come 
in  and  transfer  one  side  of  the  galvanometer-as-voltmeter  to  earth,  so  that  the 
voltage  to  earth  is  recorded.  On  the  eighth  step  the  A  group  of  relays  is  also 
brought  in  and  reconnects  the  galvanometers  as  ammeters,  R4  becomes  the  shunt 
shown  in  Fig,  6.2b,  but  for  heavier  current  R5  can  be  reconnected  via  a  39  and  a 
10  ohm  resistor  to  pick  up  l/5  of  the  voltage  across  R4. 
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6,2*4  Cycling  arrangements 

All  the  uniselectors  in  this  system  are  driven  from  a  transistor  multivibrator 
which  switches  a  power  treinsistor  via  an  intermediate  transistor.  A  master  switch 
transfers  the  current  pulses  to  temperature  scan,  leakage  scan  (described  below) 
or  power  scan.  The  temperature  scan  uniselector  is  started  by  a  press  button  and 
then  cycles  to  its  homing  position  and  stops.  The  leakage  scan  (USE),  which  has 
single  ended  wipers,  is  started  by  pressing  a  button  and  cycles  t:\rough  43 
positions;  on  its  49th  step  it  starts  USE,  which  cycles  through  49  positions  and 
on  the  50th  step  pulses  USD  so  that  both  reach  a  homing  position.  From  the  load 
scan  during  the  power  run  (already  described)  USA  and  USB  are  started  by  pressing 
a  button  and  scan  through  the  eight  positions  three  times  over  before  stopping  on 
a  homing  position.  These  uniselectors  can  be  run  continuously  by  closing  a  switch, 

6.2.5  Cesium  injection 

At  the  start  of  cesium  injection  a  button  is  pressed  which  energises  RL2, 
giving  event  marks  on  both  ultraviolet  recorders  and  on  a  recorder  used  for 
arrival  of  cesium  at  the  lower  end  of  the  heat  exchanger.  This  also  gives  a 
common  point  in  time  on  all  three  charts, 

6.2.6  Cesium  flow 

The  optical  transmission  of  the  cesium  vapour  at  the  heat  exchanger  exit  is 
monitored  by  a  simple  circuit.  The  large  window  at  this  part  is  blacked  out  and 
use  is  made  of  the  two  small  windows  in  line  with  each  other.  A  150  W  bulb  is 
placed  outside  r>ne,  while  the  other  carries  a  photocell.  An  OCP  71  phototransistor 
was  found  unacceptable  owing  to  dark  current  drift.  Next  a  CdS  cell  (NSL-33  by 
National  semiconductors  of  Canada)  was  tried.  This  is  a  specieilly  stable  type  for 
measurement,  but  its  response  to  falling  intensity  is  slow  and  it  is  being  replaced 
by  type  36-P.  The  cell  is  connected  in  a  bridge  circuit  energised  by  a  1 .5  V  cell 
and  the  output  taken  to  a  10  mV  recorder.  Ranges  of  sensitivity  can  be  provided 
by  adding  resistors  in  series  with  the  battery. 

6.2.7  Calvanometer  switching 

All  the  recorder  galvanometers  are  liable  to  be  used  for  power  measurements, 
and  panel  switches  are  fitted  to  direct  galvanometers  for  use  in  ocher  tests,  either 
with  only  a  shunt  resistor  (S3X,  34X)  or  as  a  fixed  range  voltmeter  (S1X,  S2X), 

These  switches  are  only  fitted  to  a  few  channels;  they  could  be  incorporated  in  the 
master  switch. 
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6. 2*8  Duct  Leaka/iies 


The  electrical  leakage  of  the  insulating  iaaterr'.al  of  the  duct  is  important 
and,  to  take  all  necessary  measurements  in  the  shortest  time,  a  galvanometer  ie 
used  in  a  conventional  ohmmeter  circuit  (Fig.  6.2c).  R  is  given  switched  values 
of  10,100,  3300  and  100,000  ohms,  while  the  output  leads  are  switched  via 
uniselectors  USD  and  USE.  The  readings  obtained  are  on  a  non-linear  scale,  and 
it  is  intended  to  fit  a  variable  shunt  on  the  galvanometer  to  set  FSD,  and  to 
produce  a  calibrated  scale  in  transparent  plastic  for  reading  resistance  values 
directly  off  the  chart.  When  this  system  is  in  use,  a  mask  is  used  to  black  out 
all  other  galvanometer  beams,  except  one  used  for  registering  uniselector  drive 
pulses,  which  are  picked  up  by  a  mutual  inductor  T1 .  At  fiist  an  18V  battery  was 
used,  but  the  results  did  not  correspond  with  Avo  readings  across  the  same  points; 
is  was  pointed  out  that  plasma  breakdown  was  occurring.  Accordingly,  the  battery 
voltage  was  changed  to  1.5V  emd  a  more  sensitive  galvanometer  used. 

6.2.9  Plasma  passive  electrical  characteristics 

Using  the  same  unisolector  switching,  110  V  a.c.  can  be  applied  through  a 
1  kilhom  series  resistor;  voltage  and  current  are  displayed  on  separate  galvanometers. 
This  method  gives  a  very  indirect  picture  of  the  actual  V-I  curve,  and  it  is 
intended  to  substitute  a  low-impedance  sawtooth  generator  for  the  110  V  supply. 

If  current  is  monitored  on  the  recorder  the  linear  voltage  rise  will  make  the  time 
axis  of  the  chart  correspond  to  a  voltage  axis. 

6.2.10  Electrode  resistance 

Each  of  the  36  electrodes  consists  of  a  U-shaped  loop,  and  provision  is  made 
for  connecting  both  ends  of  each  loop  to  the  rig.  When  the  master  switch  is  set 
to  ’Temperature’  contacts  are  closed  on  a  bank  of  relays  which  pass  a  cun^nt  of 
about  1A  through  all  the  loops  in  sejrvice  and  also  through  a  0.2  om  calibrating 
resistor.  Pressing  the  'Temperature  Start’  button  causes  a  uniselector  to  scan 
the  voltages  across  all  the  loops,  connecting  the  upper  electrodes  in  turn  to  one 
galvanometer  and  the  lower  electrodes  to  a  second.  For  such  low  resistances  four- 
terminal  connections  would  be  preferable,  but  tho  arrangement  gives  some  indication 
of  temperature  distribution  by  recording  resistance  changes. 

6.3  SPECTROSCOPIC  TECHNIQUES 

Spectroscopic  methods  of  temperature  sieasurement  seem  to  be  attractive  in  the 
range  1 500°K  to  2300°K.  Preliminary  investigations  into  the  spectrum  obtained 
from  cesium  suggest  that  it  may  be  possible  to  use  cesium  line  reversal  for 
measuring  temperature  in  the  duct.  Should  this  be  feasible  it  is  proposed  to 
adopt  the  optical  arrangement  used  by  Hurle,  Russo  and  Hall'  which  enables  the 
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siasultaneous  emission  and  absorption  of  the  spectral  lines  to  be  determined  using 
e.  system  of  matched  fitters  and  photomultipliers.  The  temperature  which  will  be 
measured  will  probably  be  the  average  electron  temperature  of  the  cesium  plasma 
in  the  duct. 

Spectrograms  of  cesium  from  a  standard  discharge  tube  and  a  cesium-seeded 
flame  have  already  been  obtained  and  compared.  The  temperature  of  the  flame  was 
about  1400°K  which  is  below  tiie  level  normally  studied  in  the  duct.  The  spectr>im 
showed  three  lines  of  interest:  one  in  the  green  and  two  in  the  blue  at  4555  a,nd  4595^, 
The  intensity  of  the  blue  lines  was  much  greater  than  that  of  the  green,  as  expected 
from  the  cesium  line  intensities  obtained  by  Schneider  et  al^ ,  In  the  next 
cesium-seeded  operation  of  the  IJPD  loop  it  is  proposed  to  examine  the  duct  plasma 
spectroscopically  to  assess  the  effects  of  impurities.  An  axial  view  will  be 
taken  through  the  window  in  the  head  of  the  high  temperature'  heater. 

6,4  TEMPERATURE  MEASUREMENT  DEVELOPMENT 

During  the  operation  of  the  MPD  system  errors  were  being  introduced  in  the 
-cmperature  measurements  due  to  losses  in  the  thermocouples.  The  losses  were 
mainly  attributed  to  conduction  and  radiation  in  the  thermocouple,  but  some  error 
was  probably  caused  by  thermocouple  recovery  in  the  high  velocity  gas  stream. 

This  recovery  error  was  caused  by  the  thermocouple  beir^  unable  to  follow  the  high 
frequency  temperature  changes  in  the  gas  stream,  the  effect  being  to  indicate  a 
temperature  between  the  static  and  total  temperatxire,  A  further  error  was  caused 
by  noise  introduced  in  the  thermocouple  emf  by  the  same  fast  temperature  changes, 
but  this  may  be  eliminated  by  the  use  of  suitable  low  pass  filters. 

An  Investigation  was  begxin  to  find  alternative  methods  of  temperature 
measurement,  or  to  develop  new  thermocouple  techniques.  Two  systems  based  on 
pulsed  operation  of  thennocouples  were  studied.  The  advantage  of  using  a  pulse 
system  was  that  the  thermocouple  never  reached  temperatures  where  conduction  and 
radiation  losses  become  importeuit.  This  was  brought  about  by  either  inserting 
and  withdrawing  the  thermocouple  from  the  hot  gas  at  high  rate,  or  by  operation 
with  alternate  heating  and  cooling.  The  disadvantages  of  using  this  type  of 
thermocouple  are  less  than  for  a  conventioneil  high  temperature  thermocouple,  ed though 
the  pulse  thermocouple  tends  to  read  low  with  the  rise  in  specific  beat  of  ths 
theimocouple  material  and  the  recovery  error  affects  the  thermocouple  in  the  same 
way  as  a  conventional  type. 

The  first  method  (the  ’intermittent’  pulsed  thermocouple^)  was  based  on  a 
standard  thermocouple  which  was  passed  through  a  solenoid  supplied  with  a 
rectanguleir  waveform.  A  solenoid  supplied  with  a  rectangular  waveform  vibrated 
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the  thermocoupla  so  that  it  was  imiuersad  in  the  hot  gaa  long  enough  to  reach 
the  temperatures  associated  with  high  conduction  and  radiation  losses.  The  actual 
temperature  was  then  derived  from  the  heat  transfer  between  the  gas  and  the 
thermocouple:  The  physical  aise  of  the  system  and  the  difficulty  in  determining 

the  heat  transfer  data  made  it  impractical  for  use  in  the  MPD  loop. 

The  second  system  was  made  up  from  a  chromel-alumel  thermocouple  fitted  into 
a  coaxial  tube  throxigh  which  the  coolant  flowed.  The  cooling-heating  cyc.le  was 
controlled  by  a  motorised  cam  and  microswitch.  When  the  cooling  was  stopped, 
the  thermocouple  output  rose  exponentially  (Fig.  6.3)*  Cooling  was  applied  before 
radiation  and  conduction  losses  became  significant.  The  exponential  output  was 
fed  to  a  thermocouple  amplifier  and  computer  circuit  which  predicted  the  final 
temperature.  The  limiting  factor  for  the  pulsed  thermocouple  was  the  signal-to- 
noise  ratio  of  the  thermocouple  emf;  when  reduced  to  a  certain  level  instabilities 
occur  in  the  computer  network. 

Wormser  and  Pfunter^  developed  a  practical  system  of  pulse  thermometry  upon 
which  the  present  study  was  based.  Wormser  and  Pfunter  reported  that  their  system 
was  used  to  ueasure  temperatures  up  to  2900°C  in  an  oxy-hydrogen  flsuae  with  an 
accuracy  within  1^  of  the  theoretical  temperature*  The  gas  coolant  flowed  thr:^jgh 
a  water-cooled  Jacket  directly  on  to  the  thermocouple  tip.  The  construction 
originally  proposed  was  modified  in  the  present  development  to  the  arrangement 
shown  in  Fig.  6.4. 

This  construction  tended  to  increase  the  time  constant  of  the  themao couple, 
creating  a  problem  with  the  computer  circuit  components  particularly  when  static 
gas  temperatures  were  being  measured. 

The  computer  circuit  was  composed  of  six  Solatron  Operational  amplifiers 
(Fig.  6.5)»  The  first  amplifier,A^ ,  raised  the  thermocouple  emf  to  a  value 
suitable  to  operate  the  second  amplifier, Aj ,  which  was  connected  as  a  differentiator. 
Amplifier  A^  output  had  to  have  a  large  signal-to-noise  ratio,  as  noise  fed  into 
Aj  resulted  in  instabilities  in  this  amplifier.  Aj  differentiated  the  exponential 
riso  and  added  the  input,  giving  a  steep  output  equal  to  the  actual  temperature. 


The  theoretical  evaluation  of  computer  operation  assumed  that  the  thermocouple 
emf  was  of  the  form 

-Vt 

V.J,  =  Vq  +  (1  -  e  .  (6.1) 

immediately  after  cooling  was  stopped  at  time  t  =  0. 

Tq,  and  (Fig.  6.6)  were  represented  by  V^,  and  respectively, 

where 
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Vq  =  the  cooled  thermocouple  emf  at  t  =  0, 

=  the  final  value  of  the  exponential  emf  rise, 
the  fxnal  thex'moQouple  emf ^  ^  at  t  — 

t  =  time  from  t  ”  0,  and 

r  =  the  thermocouple  time  constant,  defined  as  the  time  taken  for  the 
thermocouple  exponential  emf  rise  to  reach  0.632  V  , 

A 

The  transient  response  cf  the  amplifier  Aj  resulted  in  an  output 
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Substitution  of  equation  f6.l)  into (5, 2)  gives 
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Solving: 
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It  was  arranged  that: 


(6.5) 


and  C<R4  =  T  . (6.6) 

Therefore 

*  \  »'  '  . (6-7) 

''out  =  ''0  ^  \  . (^-8) 

This  was  the  output  of  amplifier  A3  and  was  equal  to  the  thermocouple  output  at 
t  =  00 

"■t/r 

Equation  (6,1 )  stated  \  "  ^0  ^  ®  ^ 

Let  t  =  w,  then  V  =  V  +  V  (l  -  o’”) 

X  U  A 

4  '  ''o  * 
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This  signal,  in  terms  of  voltages,  was  fed  into  the  input  of  the  third 
etmplifier  A3.  A  large  capacitor,  C? ,  wes  placed  across  the  input  and  charged  up 
by  amplifier  A2  to  a  Vv  ge  (V^  +  tlie  end  of  the  heating  cycle  the 

circuit  between  asiplifiors  A^  and  Aj  was  broken  by  RI^  /2  with  this  voltage  still 
across  Cj.  This  voltage  was  retained  at  the  input  of  amplifier  A3  by  having  a 
large  resistance,  R5 ,  in  series  with  the  first  grid  of  the  amplifier  so  making  the 
discharge  time  of  C2  large, 

Dick*  Sind  McKensie®  stated  that  r  was  dependent  on  the  physical  size  of  the 
thermocouple,  the  velocity  or  flow  rate  of  the  hot  gas,  and  the  magnitude  of  the 
step  increase  a  thermocouple  was  responding  to.  The  effect  of  R4C3  being  larger 
or  smaLler  than  t,  equation  6,6,  was  to  give  an  output  at  amplifier  A2  larger  or 
smaller  than  +  V^.  This  error  tended  towards  with  time  (Fig,  6.7)» 

To  enable  the  R4C3  value  to  be  chsinged  and  the  R4  =  R3  relationship,  equation  6.5, 
to  be  maintained,  two  ’ganged*  variable  potentiometers  were  used.  To  obteln  the 
optimum  conditions  automatically,  the  output  of  amplifier  A2  was  differentiated  by 
amplifier  A4 .  The  second  differential  of  the  thermocouple  emf  gave  the  rate  of 
change  of  the  slope  and  so  the  change  in  r, 

dV, 

For  amplifier  A4 ,  ^  CjRg 


where  given  by  equation  6*4  and  R4  =  R3 

''out  =  i  [''0  ^  ''a  -  ''a 


V  =  V  e  '  (1  - 

out  T  A  '  T 


-V. 


C.Rt 


When  CiR4  >  r,  V  .  >  V, 
*  '  out  A 


0,H.  <  r, 


C.».  =  r.  =  0 


The  output  of  amplifier  A4  was  fed  to  a  servo-amplifier  unit  Ag  which  drove 
the  two  variable  potentiometers  R3  smd  R4  to  a  position  where  R3  =  R4  and  CiR4=t. 
Under  these  conditions  the  second  differential  of  the  thermocouple  emf  was  zero. 

By  a  system  of  relays  (Fig,  6,8)  this  operation  was  timed  to  take  place  during  the 
cooling  cycle,  as  movement  of  Rj  and  R4  during  the  heating  cycle  introduced 
instabilities.  The  servo-amplifier  unit  Ag  introduced  a  great  deal  of  noise  into 
the  system  and  so  an  integrator  A5  was  connected  between  amplifiers  A4  and  Ag , 
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Amplifier  A,  was  suitably  switched  to  operate  during  the  heating  cycle. 

Up  xo  the  present  time,  only  amplifiers  ,  Aj  and  A,  have  been  successfully 
used  «s  it  has  not  been  possible  to  eliminate  all  the  instabilities  in  amplifiers 
A4 ,  Aj  and  Ag .  The  circuit  used  was  sufficient  to  measure  temperatures  when  t 
did  not  change  rapidly  6u^d  '^■his  was  done,  although  some  inaccuracy  probably  occurrea. 

Two  sizes  of  thermocouple  were  investigated  (Fig.  6.9);  the  larger  is  0,062  in. 
diameter  mineral  insulated  wire  type  and  the  smaller  a  C^OA-  in.  diameter  mineral 
insulated  wire  type. 

Under  static  gas  conditions,  T  varied  considerably  with  temperature,  Fig. 6. 10. 
The  value  at  r  was  large  for  the  0.062  in.  thermocouple,  i.e.  “  0.110  sec  at 
300^C,  This  large  value  of  t  required  large  C1R4  values  for  amplifier  Aj .  The 
type  of  amplifier  used  was  not  suited  for  very  large  component  values  such  as  these. 
With  a  value  of  R4  >5  Mohm,  without  noise,  the  amplifier  was  unstable,  6md  even 
with  this  value  a  large  was  required  i.e.  22^f  which  created  overloading  problems 
in  the  amplifiers.  The  system  was  calibrated  up  to  800°C,  using  manual  adjustment 
for  T  =  C, R4  conditions.  The  repeatability  was  good,  but  the  instabilities  in 
amplifier  created  oscillation  about  the  output,  making  accurate  calibration 

difficult. 

The  use  of  the  0.04.  in.  thermocouple  reduced  t  by  a  factor  of  two  at  300*^C, 
and  at  800°C  t  was  approximately  20  sec. 

From  Fig.  6.10  it  was  seen  that  r  fell  with  increasing  temperature  and  from 
this  trend  the  required  C<R4  should  be  of  realistic  values  at  temperatures  above 
lOOO^C. 

Only  preliminary  tests  have  been  carried  out  with  flowing  gas  to  deteriiilne  the 
magnitude  of  t.  Using  an  ordinary  gas  torch,  t  was  50  sec  at  500°C  and  12  sec  at 
1100°C  (Fig.  6.11)  tor  the  0.062  in,  thermocouple.  Again  it  was  assumed  that  as 
the  temperature  rose  above  lOOO^C  the  value  of  r  would  reduce  even  further.  This 
reduction  of  r  and  C1R4  would  give  Increased  stability  and  reduce  the  tendency  for 
amplifier  Aj  to  overload, 

A  more  stable  system  should  be  developed,  perhaps  by  using  servo-differentiators 
is  they  were  found  to  have  any  advantages.  It  is  proposed  to  extend  the  range  of 
measurement  to  at  least  2000°C  under  flow  conditions,  and  also  with  static  gas  if 
feasible.  The  high  value  of  t  was  still  troublesome  and  therefore  a  new  design 
of  coolant  tube,  to  reduce  t,  is  to  be  considered;  with  this  reduction  in  r 
achieved,  the  instabilities  in  amplifier  A2  would  probably  be  eliminated.  The 
feasibility  of  using  one  computer  system  for  several  thermocouples  is  also  to  be 
studied. 
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(a)  qALVANOMCTCR  AS  VOLTMETER. 


(b)  GiALVANOMETER  AS  AMMETER 


(c)  GALVANOMETER  AS  OHMMETER. 
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THERMOCOUPLE  RESPONSE  WHEN  COOLED 
AUTOMATICALLY  AT  SET  TIME  INTERVALS 
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FITTIJRE  PROGRAMME 

by 

B.C.  Lindley 

As  intimated  in  the  previous  Technical  Summary  Report  (IRD  64**36),  a 
continuing  two-year  programme  of  research  is  envisaged,  following  the  successful  1 

technological  development  which  led  to  the  first  cesium-seeded  operation  in 
Maroh  1964*  Much  of  the  value  of  the  data  obtainable  in  subsequent  operations 
depends  upon  ability  to  obtain  the  f]  w  parameters  with  reliable  aoouraoy  and  to 
eliminate  the  leakage  currents  due  to  electrically-conducting  layers  on  the  MFD 
duct  walls.  Convenience  of  operation  of  the  facility  will  depend  on  suooess  in  | 

developing  methods  of  separating  cesium  from  helium,  a  problem  which  is  important  I 

in  a  number  of  practical  concepts  of  MPD  generator. 

Details  of  the  range  of  parameters  over  which  data  on  plasma  behaviour  and 
MFD  power  extraction  are  to  be  obtained  were  given  in  Section  1  .2  of  IRD  64-36. 

A  somewhat  greater  emphasis  will  be  placed  on  electrode  processes  since  it  is 
believed  that  these  represent  an  ultimate  limitation  on  specific  power:  current 
densities  measured  to  dile  Have  been  up  to  a  few  hundred  and  in  large-scale 

systems  two  orders  of  magnitude  Increase  may  be  necessary.  | 

I 

I 

I 
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CHAPTKR  8 


KICROWAVB  arUDIES  OF  NON-SQUILIBRIUM  PLASMA 

by 

J.H.  H&rrls  asd  D.  Balfour 

8.1  SBNSRAL 

The  microwave  cavity  experimert  was  originally  suggested  as  a  means  of 
studying  photolonlsatlon  and  recombination  in  cesium  and  oeslum-helitus  mixtures 
in  relation  to  the  possibility  of  increasing  the  efficiency  of  MPD  generators  by 
using  ultra  violet  radiation  to  induce  non*^quillbrl\im  ionisation.  This  reskains 
the  long-term  objeotlve  but  in  testing  this  experiment,  in  understanding  the 
behaviour  of  the  microwave  bridge  and  in  calibrating  the  response  of  the  bzddge  it 
was  advantageous  to  use  eleotrioal  discharge  pulses  to  create  the  ionisation. 

Using  a  discharge  tube  it  was  possible  to  out  off  the  electric  field  in  a  few 
mioroseoonds  and  then  study  the  loss  processes  occurring  in  the  afterglow  of  the 
discharge. 

The  information  which  this  experiment  can  yield  on  the  decay  prooesses  in 
oeslum-heli\uD  mixtures  is  of  very  great  importance  to  MPD  generation  in  that  the  ! 

degree  of  ionisation,  and  henoe  the  eleotrioal  conductivity  attainable  in  an  MHD 
generator  under  non  equilibrium  conditions,  is  limited  by  the  rate  of  loss  of 
ionisation.  The  loss  prooesses  include  diffusion,  recombination  and  attachment.  i 

The  diffusion  luss  may  be  small  in  a  generator,  and,  while  attachment  to  cesium  | 

fmd  helium  may  be  negligible,  attachment  to  impurity  atoois  may  be  slgnlf leant. 

i 

A  small  quantity  of  oxygen  for  example  could  make  attachment  the  dominant  loss 
mechanism.  Most  theoretical  papers  on  generator  performanoe  have  considered 
recombination  as  the  only  loss  process  and  in  particular  have  used  values  of 
recombination  coefficient  of  •'*  10”^*  omVion  seo  appropriate  to  a  three-body  process. 
However,  at  lower  temperature,  dissociative  recombination  nay  be  more  important  and 
the  value  of  10**  cmVlon  sec  measured  by  BLondi  and  Brown' was  ascribed  to  this 
process.  The  uncertainty  of  the  recombination  mechanism  and  the  very  sensitive 
dependence  of  three-body  recombination  on  temperature  call  for  accurate  measirements 
of  the  reoGoblnatlon  coefficient. 

AH  of  these  loss  prooesses  can  be  studied  in  the  existing  experiment.  So 
far  the  pressure  and  degree  of  ionisation  have  shown  diffusion  to  bo  the  essential 
loss  process;  measurements  have  been  made  of  the  mobility  of  cesium  ions  in  cesium 

I 

vapour.  The  ionisation  and  pressure  can  be  adjusted  to  enable  reoombination  to  I 
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be  more  important  and  the  measurement  of  recombination  of  cesium  will  bo  the  next 
stage  in  the  experimental  progrEunmo..  A  further  stage  will  bo  the  introduction  of 
impurities  to  oteasxire  attuohment  probabilities* 

8.2  MICRCWrAVK  BRIDSK 

The  previous  Technical  Summary  Reports  63-1  and  6if-36)  describe  how  a 
steel  cavity  was  designed  and  constructed  to  operate  in  the  TM020  i&ode  at  a 
resonant  fre<pjenoy  of  approximately  9000  mc/sec  with  a  quarts  tube  Inserted  along 
its  axis*  The  cavity  was  incorporated  in  a  microwave  bridge  circuit  (Pig.  8.1 ) 
or  frequency  discriminator,  whore  a  change  in  resonant  frequency  of  the  cavity 
results  in  an  out-of -balance  signal* 

Tests  on  this  bridge  were  made  usln/':  the  ionisation  produced  in  an  electrical 
discharge  in  air  between  tungsten  electrodes  seeled  in  the  ends  of  the  plasma  tube. 
The  first  type  of  discharge  to  be  used  was  produced  by  half-wave  rectifisd 
alternating  high  voltage  supply*  This  gave  results  of  the  type  illustrated  in 
Fig*  8*2  where  the  top  tz^ce  shows  part  of  one  cycle  oi  the  current  flowing  through 
the  discharge  tube,  azui  the  bottsa  truce  gives  the  out-of -balance  signal.  A  square- 
wave  modulation  of  lOCCc/sec  was  applied  to  the  klystron  signal  in  this  early  set 
of  measurements  and  both  traces  are  displayed  with  time  plotted  horizontally.  Tb* 
interesting  feature  of  Fig.  8*2  is  that  over  tL<  region  whezc  the  c>  '-.-ent  is 
inozH»asing,  the  bridge  signal  first  grows  z^pidly,  passes  through  a  maximum  and  then 
slowly  decreases* 

8*3  BRIDGE  THKORY  AMD  CALIBRATION 

A  theory  of  the  action  of  the  bridge  (Fig*  8*1 )  was  developed  to  explain  these 
features,  where  the  bridge  has  been  adjusted  initially  so  that  in  sirm  A  the  signal 
C  reflected  ft*om  the  cavity  leads  the  signal  P  reflected  from  the  jshort-circuit 
plunger  by  90°  but  in  arm  B  the  signal  C  lags  90°  behind  the  signal  P*  This  is 
Illustrated  in  Fig*  8*3* 

Bridge  out  of  balance  signal  -  B^  **  A^ ,  since  square  law  detectors  are  used, 

»  icPC'sin  t4> 

where  C*  is  the  slgzud  reflected  fr^  the  cavity  and  is  the  shift  in  phase  of 
the  reflected  signal  when  free  electrons  are  present  in  the  plasma  tube. 

If  7  is  the  voltage  amplitude  falling  on  the  cavity  and  S  is  the  VuffR 
of  the  cavity,  then 

Bridge  signal  =  itPV  (^)  sin  .  (8.1 ) 

The  reactance,  X,  of  the  cavity  is  given  by 

X  =  2Q^R  ^  .  (8.2) 
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where  E  ia  the  cavity  reaistance, 

ia  the  unloe  led  aeleotivity  of  the  cavity, 

and  Af  ia  the  ahl'  -  of  frequency  from  the  resonant  froq’  ^ncy  f^. 

With  a  new  quartz  plasma  tube  inserted  in  the  cavity, the  precision  attenuator  and 
alotted-guide  standing  wave  indicator  were  used  to  obtain  values  of  S  over  a  range 
of  frequencies  with  the  results  shown  in  Fig*  8*4*  The  graph  yields  the  following  values: 

VSWR  at  resonance  =  0*50 

loaded  aeleotivity  (Qj^)  =  5700 

unloaded  selectivity  (Q^)  =  8000  and 

coupling  factor  (fi)  =  0*40 

The  value  of  S  at  resonanoe  specified  the  impedance  contour  to  be  used  on 
a  Smith  Chart.  A  value  of  Af  was  chosen  and  substituted  in  equation  (8*2)  to 
obtain  the  coz*respondlng  X  value.  The  intersection  of  this  value  of  X  with  the 
impedance  contour  gives  the  corresponding  values  of  both  A^  S  to  be  substituted 
in  equation  (8*1 )  to  calculate  the  bridge  signal.  This  procedure  is  Illustrated 
in  fig.  8.5.  Results  were  compared  with  two  independent  sets  of  experimental 
results.  One  method  involved  keeping  the  klystron  frequency  constant  and  varying 
the  resonant  frequency  of  the  cavity  by  varying  its  diameter  (throu^  change  of 
teaperature  ever  a  range  of  40°C).  This  gave  the  results  shown  in  Fig.  8.6.  In 
the  other  method  the  resonant  frequency  of  the  oavity  was  kept  constant  and  the 
klystron  frequency  was  varied;  this  gave  an  identical  set  of  results.  Fig.  8.6 
shews  the  agreement  between  the  theoretical  and  experimental  results.  The 
variation  of  out-of ‘■balance  signal  with  frequency  shift  is  called  the  calibration 
curve. 

In  some  experimental  oases,  the  calibration  curve  was  observed  to  be 
asymmetrical.  This  may  possibly  have  been  due  to  inaoourate  setting  of  the  short 
circuit  plunger,  and  to  check  on  this,  the  effect  on  the  theoretical  curve  of  an 
inaccuracy  of  setting  of  1mm  was  evaluated  by  modifying  equation  (8.1 )  to  beoome: 

Bridge  signal  =  4PV  .....  (8.3) 

where  a  is  the  phase  change  introduced  by  a  plunger  movement  of  1  mm.  The 
resulting  calibration  curve  is  shown  in  Fig.  8.7. 

The  decay  of  ionization  on  suddenly  removing  the  voltage  applied  to  an 
electrical  discharge  in  air  was  recorded  using  different  sections  of  the 
calibration  curve.  The  bott<»i  trace  of  Fig.  8.8  shows  a  plasma  decay  taking  the 
cavity  over  the  range  A  to  D,  Fig.  8.6b  from  B  to  3  and  Fig.  8.8o  from  C  to  F 
(for  the  significance  of  the  letters.  A,  B,  C,  D,  B  and  F,  see  Fig*  8*9).  In  each 
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case  the  top  '■•race  shows  the  current  through  the  discharge.  Pig.  8,10  shows  the 
electron  concentration  decay  our7e8  calculated  from  Pigs.  8.8a,  b,o  and  since 
these  exhibit  the  same  form  the  interpretation  of  tho  calibration  curve  was 
considered  valid. 

8.4  IDHfTIFICATION  OF  TMni>n  MODS 

It  was  observed  experimentally  that  on  inserting  the  plasma  tube  through 
the  cavity  tho  frequency  of  tho  TM020  mode  was  shifted  by  800  Mc/s  and  that  with 
the  tube  in  position  this  resonance  was  within  60  Mc/s  of  another  mode  (see  IBD 
64-36  Pig*  16.1).  Initially  there  was  difficulty  in  identifying  the  TM020 
Inconolusivo  attoicpts  to  distinguish  the  TM020  mode  were  made  firstly  by  heating 
the  oavlty  and  secondly  by  Introducing  shims  to  increase  the  cavity  length,  but 
the  resonant  frequencies  of  both  modes  were  Independent  of  the  length  of  the  cavity. 
However,  by  inserting  a  steel  wire  along  the  axis  of  the  cavity  one  mode  exhibited 
a  continuous  decrease  in  Q  value  as  the  wire  was  Introduced  and  finally  disappeared 
completely,  whereas  the  other  mode  was  unaffected.  The  affected  mode  was  the  TMoac 
with  its  strong  axial  electrio  field. 

8.5  BiaCTRON  CONCBNTHATIOK  MEASURafija*lTS  IN  ARGON  APTfiRSLOW 

Having  established  the  working  charactsristlcs  of  the  bridge  and  its  ability 
to  follow  a  rapid  decay  in  electron  concentration,  attention  was  turned  towaird 
decay  conditions  in  argon  to  obtain  a  comparison  with  results  reported  by  other 
workers.  The  air  in  the  discharge  tube  was  replaced  with  commercial  cylinder 
argon  and  the  decay  of  ionisation  on  suddenly  removing  the  applied  voltage  was 
rsoorded  over  a  ran^o  cf  argon  pressures. 

Pig.  8  •11  shews  a  typical  oscilloscope  trace,  where  the  argon  pressure  was 
2.7  torr:  the  graticule  squares  each  represent  20  fiseo  on  the  horizontal  time 
scale.  The  top  trace  shows  the  bridge  signal,  line  1  gives  the  level  of  the 
upper  turning  point  of  the  calibration  curve  and  line  2  gives  the  lower  turning 
point;  the  balance  point  of  the  bridge  in  this  trace  was  on  the  lower  turning 
point.  Tho  variation  of  electron  concentration  with  time  calculated  from  Pig. 8. 11 
is  shown  in  Pig.  8.12  (the  calculation  made  to  obtain  electron  concentration  from 
frequency  shift  assumes  that  the  collision  frequency  is  small  compared  with  the 
radian  frequency  of  the  applied  signal). 

Heotron  loss  processes  in  the  decaying  plasma  are  ambipolar  diffusion, 
ion-eleotrcn  recombination  euid  electron  attachment,  and  the  loss  rate  may  be 
written  as 

R  D  7^  n  -  an  *  -  hi;  n 

a  e  e  0  e 
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dn 

_ e 

dt 


(8.4) 


where  is  the  electron  concentration. 

Is  the  anblpolar  diffusion  coefficient, 

a  is  the  recombination  ooefficierit, 

h  is  the  attachment  coefficient,  and 
is  the  (dleotron  collision  ft*aquency, 

A  study  of  equation  (8^4)  shows  that  the  diffusion  and  attachment  prootsres 
yield  an  exponential  decay  of  electron  concentration  with  time,  but  the  recombination 
proeeas  leads  to  a  linear  relationship  between  the  Inverse  of  eleetron  concentration 
and  time# 


In  fig*  d*13  the  logarithm  of  frequency  shift  is  plotted  «g;ainet  time«  8iaee 
frequenoy  shift  is  proportional  to  sleotron  aenoentratlon  tlie  straight  line  ton  of 
these  ourves  indicates  that^as  an  electron  loss  prooess^reoMbination  is  negligible 
at  the  gas  pressure  oxhSl  elootron  concentration  used  in  the  present  experiments* 
furthezmore,  the  rate  of  loss  of  eleotrons  varies  approximately  as  the  Inverse  of 
the  gas  pressure  and  this  suggests  that  the  loss  process  is  essentially  diffusloni 
the  pressure  dependence  of  an  attaohment  process  would  be  either  linear  in  ^the  oase 
of  a  two«d>ody  meohaniun  or  quadratic  in  a  thjree-body  mechanism*  A  theoretical 
treatment  of  ambipolar  diffusion  leads  to  tho  following  expressions  for  the  ambipolar 
diffusion  ooeffioient  and  the  reduced  mobilily  (reduced  to  RCP) 

Da  ■  ^  ^  (^0**  firet  order  diffusion)  •••••  (8*5) 

■  ^  ^  om*/volt  sec  . (8.6) 

1  i 

T 

« 

where  r  is  the  radius  of  tho  disoharge  tidbe, 
p  is  the  gas  pressure  in  torr, 

A  is  the  characteristic  diffusion  length  of  the  plasma  enolosure, 
r  is  the  tine  constant  of  tho  electron  concentration  decay,  and 
is  tho  gas  temperature  and  the  eleotron  temperature. 

In  tho  results  qucted  hero,  the  eleotron  temperature  was  assumed  to  be  equal 
to  the  gae  temperatuire.  When  the  applied  voltage  is  removed  the  eleotron  temperature 
ie  likely  to  be  much  larger  than  the  gas  temperature  but  it  has  been  estimated  that 
thermal  eq\>ilibrlum  will  be  approached  in  a  period  of  approxiaately  30  gaeo*  The 
sero  of  time  in  Figs.  8.11,  8.12  and  8.13  oorrespoaid  to  at  least  30  n»oo  after  the 
removal  of  the  ionising  voltage. 


e.5 


lilllWimiillillliiWIWIlB'liillllWIIIWlWHI'l'*""'' . . . . .  '■““•""“"‘"“m'HMHilk.wiili*,  i-'4i|i,|jt|(i||ii(i()'|(|il|ii|i|iliHIW!«2le, 


Bioxidi  and  Brovn  haT«  shown  that  for  ambipoleu?  diffusion  under  oonditioixa  of 
thermal  equilibrium  *  oonatani,  where  is  the  pressure  reduced  to  normal 

teaperatura.  Values  of  and  p^  for  argon  orer  the  range  of  pressures  studied 
were  calculated  from  the  ourres  of  Fig*  8*13  and  are  gi-yen  in  the  following  table. 


p,  tOXT 

0.1  1*5  2*7  6*0  8.2  10 

omVvolt-seo 

p  D  ,  torr  cffi^/seo 

0 

7.9  9A  4*9  6.4  7.5  8*3 

350  420  220  260  330  370 

The  values  of  p  and  D  are  suffioient3.y  constant  to  bo  consistent  with  electron 
loss  duo  tu  rigblpolar  diffusion  as  deduced  independently'  from  the  variation  of  the  loss 
rate  with  pressui'e*  The  values  of  shown  in  the  above  table  stb  higher  by  a 
faoto*  of  at  least  2  than  thr  values  reported  by  Biondi  and  Chanin^*  Howeveri  the 
impurity  in  the  argon  in  these  experiments  may  have  been  appreciable*  Higher  values 
of  elaotron  ooncentration  are  required  before  recombination  beoomes  the  dominant 
prooeas* 

8*6  BLBCTROy  MUaiTRTlfgNTa  IK  CESIUM  AFTHR&LCBT 

The  decay  prooesa  in  oesiun  vapour  was  studied  in  the  same  way  as  for  argon. 

The  eleotrioal  discharge  tube  was  connected  to  a  cesiufl»cap8ule  container  and  bath, 
pumped  down  to  a  very  low  pressure,  sealed  *f  and  enclosed  in  o^ens*  The  capsule 
was  broken;  the  plasma  tube,  connecting  line  emd  cavity  were  raised  to  a  constant 
temperature  in  the  region  of  TOO^X^and  the  oesium  vapour  pressure  was  controlled 
and  e8tima»ed  using  the  temperature  of  the  cesium  bath*  High  voltage  piilses  were 
applied  to  the  cesium  vapoxir  and  out  off  sharply  using  the  method  described  in  the 
previous  report  (IHC  6i4.“36)* 

The  ixiitial  experiments  were  aimed  at  studying  the  effect  of  the  power  of 
the  probing  signal  on  the  decay  of  electron  concentration,  since  the  microwave  signal 
might  tend  to  heat  the  elec  crons «  The  effect  of  varying  the  length  of  the  ionisation 
pulse  so  that  uniform  breakdown  might  be  produced  in  the  discharge  tube  was  also 
studied* 

8*6.1  liTfeot  of  varying  microwave  power 

Decay  curves  were  recorded  first  with  the  full  klystron  power  (appiroximately 
0*15  watt)  aupplici  to  the  bridge,  then  with  6  db  attenuation  between  klystron  &nd 
bridge, and  finally  with  13  db  inserted.  The  trace  with  $  db  inserted  is  showa  in 
I'ig*  8*14  together  with  startir^g  level,  top  bottom  letels  of  the  calibration 
curve  and  the  disohart^s  ourrent  trace* 
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The  graphs  of  the  logarithm  of  resonant  frequency  shift  of  tlie  cavity  against 
time  are  given  in  Fig*  8*1 5*  The  relative  displacements  of  the  parallel  lines  are 
due  only  to  different  times  of  triggering  and  the  straight  line  form  shoirs  that  the 
process  responsible  for  decay  of  electron  concentration  is  ambipolar  diffusion. 

The  graphs  also  show  a  diffusion  rate  independent  of  power  level. 

The  effect  of  introducing  large  microwave  power  into  the  discharge  would  be 
to  raise  the  electron  temperature  and  since 

\  ■  . (8.7) 

1  e 
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the  ambipolar  diffusion  coefficient  would  be  correspondingly  affected*  The 
measurements  indicate  that  for  the  values  of  microwave  power  used  in  tne  experiments 
the  heating  of  the  electrons  was  not  significant.  However,  under  ccudltlons  where 
recombinatiion  is  the  dominant  loss  process  the  dependence  on  electron  temperature 
may  be  much  more  sensitive.  Under  these  ciroumstanoes  a  further  check  on  the  effect 
of  microwave  pov'er  will  be  necessary. 

8.6.2  Effect  of  varying  the  length  of  the  lonlaatlon  pulse 

Various  authors  have  reported  that  the  rate  of  loss  of  electrons  trcm  afterglows 
is  dependent  on  the  duration  of  t'  a  ionising  pulse.  In  certain  oases  where  the 
ionisation  pxilse  ;..s  short  the  distribution  of  ionisation  in  the  discharge  tube  jjs 
inhomogeneous;  Oskam  and  Mittelstadt^,  for  example,  have  reported  asynmetrloal  light 
distribution  for  ionisation  pulaes  shorter  than  15  lisoo*  Where  the  ionisation  is 
limited  to  a  filament  within  the  discharge  tube  the  assiimption  that  the  ionisation 
flUt  the  discharge  tube  will  lead  to  overestimates  of  diffusion  coefficients  and  of 
recombination  coefficients.  By  increasing  the  width  of  ^  ionisation  pulse  a  more 
uniform  distribution  of  the  ionisation  can  be  azhieved  and  more  reliable  values  of 
diffusion  cmd  rsccmbinatlon  ooefficlents  can  be  deduced. 

An  experiment  was  carried  out  to  make  certain  that  the  results  deduced  from  the 
afterglow  measurements  were  not  affected  by  the  duration  of  the  ionisation  pulse* 
Cesl'mi  vapo’or  was  contained  in  the  discharge  tube  at  600°K  and  0.6  torr  as  determined 
by  the  tomperaturo  cf  the  vapour  bath.  The  duration  of  the  ionisation  pulse  was 
varied  between  100  and  800  ijsbo  and  for  each  case  the  decay  in  electron  concentration 
after  cutting  off  the  ionising  pulse  was  monitored  by  reoox^'lng  the  out-of -balance 
signal  of  the  miorowave  bridge. 

i.s  in  the  earlier  experiments  the  electron  loss  prooess  was  ^minated  by 
ambipolar  diffusion  and  since  the  frequency  shift  of  the  cavity  due  to  ionisation  in 
the  discharge  tube  is  proportloxuil  to  electron  concbintration  a  linear  dependence  was 
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found  between  the  logarithm  of  tlia  frequency  shift  and  time  (Fig,  8.1 6),  Points 
are  shown  for  the  two  extreme  values  of  ionisation  pulse  width  and,  within 
experimental  error,  all  points  can  be  related  to  the  same  straight  line. 


The  conclusion  from  this  experimental  result  is  that  ionization  pulses  of 
only  100  /isec  8u:'e  sufficient  to  produce  homogeneous  breakdown. 

8.7  CESIUM  ION  MOBILITY 

The  mobility  of  the  cesium  atoms  can  be  deduced  from  the  measured  diffusion 
coefficient.  The  Einstein  relationship  between  the  mobility  and  diffusion 
coefficients  of  a  particle  can  be  written  as 


D 


e 

kT 


(8.8) 


where  ft  is  the  nobility  of  the  particle 

D  is  the  diffusion  coefficient  «:>.f  the  particle 
T  is  the  temperature  of  the  pcrticles 
e  is  the  electronic  charge  and 
k  is  Boltzmann's  constant 

For  ions: 

r  -  ^  . 

+  + 

» 

and  the  ion  diffusion  coefficient  oan  be  expressed  in  terms  of  the  axspxbolar 
diffusion  coefficient,  D  ,  by  the  equation 


•  0+  (1  *  r  )  , 
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where  T  and  T  are  the  electron  and  ion  temperatures. 
Hence  n 


e 

kT 


ft 
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In  practice,  mobilities  aro  reduced  to  normal  temperature  and  pressxire  so 
that  the  reduced  mobility  1®  given  by 


(8.12) 


where  p  is  the  gas  pressure  in  tcrr  and  T 


g 


is  the  gas  temperature. 
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To  deduce  the  aaibipolar  diffusion  coefficient  the  diffusion  equation  must  be 
solved  for  the  cylindrical  geometry  of  the  discharge  tube. 

^  -  dlv  (ot) 

=  D  V^n  . (8.13) 

ft 

If  the  electron  concentration  decays  with  a  time  constant  r,  then  n  »  n^e  ^ 
whero  n^  is  the  initial  electron  concentration,  and  the  diffusion  equation  can  be 
written  as 

y*  Ho  +  =  0  . (8.14) 

A 

For  a  cylinder  of  radius  r^  and  height  h  a  solution  of  this  equation  is 

n(r,s,t)  =  Z  S  G  J  (a,,  r)  cos  ®xpr~ 

i»1  j=i  1  n  j 


where  is  the  i^^  root  of  the  Bessel  function  J  .  Assuming  that  the  first 

10  o 

diffusion  mode  dominates  the  loss  process,  then  i  s  j  s  o  and 


I 


For  the  discharge  used  in  the  experiment  h  » 
temperature  is  equal  to  the  gas  temperature  T 
written  as 


so  that,  assuming  that  the  ion 
,  the  reduced  mobility  can  be 


|-r  ^2^405^  f  ®“*Aolt  sec 


T 

1  +  ;;r 
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Frcm  the  slopes  of  the  log  Af  vt  plots,  walues  of  l/^-  were  deduced  for  the 
I  experiments  described  earlier  in  which  the  microwave  power  was  varied  and  in  which 
the  width  of  the  ionisation  p\ilse  was  changed.  The  calculated  values  of  deduced 
from  these  results  are  given  below. 


1  p,  torr 

T  °K 

7»  sec’^ 

cmV^olt  sec 

Experiment  1 

1.1 

630 

460  X  10’ 

0.26 

Experiment  2 

0.6 

600 

4.84x  10’ 

0.17 
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Bftrller  aeasiirementa  of  the  mobility  of  ceBium  ions  in  cesium  vapour  are  as 
follows. 


ion 

Method 

Dandurand  and  Holt^ 

Ca2^ 

0,065 

Diffusion  :  microwave  cavity 

Chen  and  Raether^ 

Cs"^ 

0.4(revised  to  0.2) 

Diffusion  ;  microwave  transmission 
(assumes  Blancs  Law  for  mobility 
in  mixtures) 

Chanin  and  Steen* 

Csj'^ 

0.21 

Time  of  flight 

Cs"^ 

0.075 

The  '^ues  of  reduced  mobility  determined  in  the  present  experiment  are  close 
to  the  revised  value  of  Chen  and  Raether  and  the  faster  ion  detected  by  Chanin  and 
Steen.  There  is  no  experimental  evidence  whereby  the  ions  dominating  the  diffusion 
loss  In  the  present  experiment  may  be  identified  but  the  values  measured  are  close 
to  the  theoretical  value  of  0.2  deduced  by  Chanin  and  Steen  from  the  Langevin 
classical  theory  of  mobility  for  Csj^. 


It  is  interesting  to  note  that  Oakam  and  Mittelstadt^  measured  veJ.ues  of 
mobility  for  rare  gas  ions  in  the  parent  gas  which  increased  with  increasing 
pressure  before  reaching  a  constant  value  at  pressures  above  a  few  torr;  this  low 
pressure  behaviour  was  explained  in  teruu  of  diffusion  cooling.  A  similar 
behaviour  might  be  expected  in  cesium  vapour  and  so  it  is  proposed  to  measure 
as  a  function  of  pressure  in  the  near  future. 

Brown^  has  shown  that  in  a  d.c.  electric  field  the  mobility  of  an  ion  oan 
bo  written  as 


m 


m. 


(8.16) 


where  is  the  mass  of  the  ion  and  is  the  frequency  of  collision  for  momentum 
transfer. 

(1  -  ooo6)  . (8.17) 

whore  v  is  the  collision  frequency  and  1  -  cos  $  is  averaged  over  all  values  of 
6f  the  angle  through  which  a  particle  is  scattered  with  respect  to  the  eleotrlo 
field  direction. 


The  total  collision  cross  section  for  the  ion  may  be  expressed  in  terms  of 
the  collision  frequency  as 

» 
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trm 
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An  accurate  asseasrent  of  1  >  oob6  la  difficult  to  determine  but  If  the 
value  1  be  uaed  then  a  lower  limit  of  500  x  10"^^  la  obtained  for  the  total 

croaa  aeotlon  for  acatterlng  of  oealum  Iona  by  cealum  atoma. 

8.8  DIFFUSION.  BSCOMBINATION  AMD  ATTACHMSNT 

In  all  of  the  experlmenta  reported,  the  loaa  of  electrona  In  the  afterglow 
haa  been  due  to  amblpolar  dlffualon*  To  meaaure  recombination  coefflclenta  It 
will  be  neceaaary  to  modify  condltlona  ao  that  the  rate  of  volume  recombination 
becomea  greater  than  the  rate  of  dlffualon  loaa.  Gray  and  Kez*r*,  ualng  a  THoio 
microwave  cavity,  have  atudled  the  decay  of  lonlaatlon  in  which  dlffualon  and 
volume  recombination  are  th  algnlflcant.  Th^  define  a  parameter  which  la  the 
ratio  of  the  initial  volume  recombination  rate  In  the  abaenoe  of  dlffualon  to  the 
rate  of  dlffualon  in  the  fundamental  mode,  and  for  reliable  eatimatea  of  recombination 
to  be  deduced  from  ^  vt  plota  they  have  ahown  that 

/?  >  15 

where  ^  = 

and  A  la  the  character  la  tic  dlffualon  length  (A  «  2 .405  for  the  dlacharge  tube  uaed 
in  the  preaent  experiment). 


Tbua  it  may  be  poaalble  In  the  preaent  apparatua  to  meaaure  loaa  d\ie  to 
volxime  recombination  by  Increaalng  the  electron  concentration  and  by  working  at 
higher  preaaurea  where  D^  la  amaller  (beoauae  p  D^  =  oonatant).  S<»a  Inoreaae  In 
electron  concentration  can  be  achieved  by  decreaaing  the  value  of  the  current 
limiting  realator  amd  the  oealum  vapour  preaaure  can  be  increaaed  to  about  10  torr 
In  the  present  apparatua.  The  diffusion  loaa  could  alao  be  reduced  by  applying  a 
longitudinal  magnetic  field  azid  the  dlmenalona  of  the  tube  could  be  increased  if 
necessary.  Thus,  by  choosing  suitable  conditions  It  should  be  poaalble  to  enable 
]reocmd>lnatlon  coefficients  to  be  measured  In  pire  cesium.  Subsequently  the 
recombination  processes  occurring  In  oealum-helium  mixtures  will  be  studied  and,  by 
Introducing  impurities  Into  the  plasma  tube,  it  will  alao  be  possible  to  evaluate 
attachment  coefficients. 
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The  experiicent  ia  inherently  capable  of  studying  all  the  loss  processes  that 
will  limit  the  level  of  ionisation  and  plasma  electrical  conductivity  attainable 
in  an  MPD  generator. 

8.9  PRELIMINARY  EXPERIMENTAL  STUDIES  OF  PHOTOIONIZATION 

Fig.  8.17  shows  a  plasma  tube  with  quartz  end  windows  which  was  used  to  detect 
photoionization  in  oesium  vapour  produced  by  a  merc\uy  vapour  lamp.  The  cesium 
vapour  pressure  was  determined  by  controlling  the  temperature  of  the  cesium  bath 
and  the  plasma  tube  was  heated  to  approximately  600°K.  The  variation  of  electron 
density  in  the  ces^'jm  vapour  is  shown  in  Pig.  8.18  for  a  vapour  pressure  of 
approximately  10  torr  and  with  the  lamp  mounted  so  that  an  ultra  violet  beam  was 
directed  along  the  axis  of  the  plasma  tube.  The  lamp  was  operated  from  90  c/s 
mains  supply  and  the  electron  density  reveals  the  anticipated  100  c/s  variation 
superimposed  or  a  much  lower  frequency  which  was  due  to  a  mechanical  vibration. 

This  lower  frequency  is  also  displayed  in  the  absence  of  any  illumination  from  the 
lamp. 

The  microwave  bridge  was  tuned  to  maximum  sensitivity  to  record  this 
photoionization  and  the  maximum  electron  density  detected  was  approximately  ic’  cm~^. 
The  variation  of  peak  electron  concentration  as  a  function  of  vapour  pressure  is 
shown  in  Fig.  8.19*  The  difference  between  the  results  for  heating  and  cooling  is 
due  to  a  lag  between  the  cesium  vapour  pressure  variation  and  the  bath  temperature 
variation. 

These  measurements  are  a  satisfactory  preliminary  step  in  the  study  of 
photoionization  but  for  measurements  of  photoionization  cross  section  and  for 
recombination  a  more  powerful  source  of  ultra  violet  radiation  will  bo  required. 
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MICROWAVE  BRDGE  OR  FREQUENCY  DISCRMtSWTOR 


CURRENT  THRCXiGH  DISCHARGE  TUBE  AND  CORRESPONOWC  OUT  OF  BALANCE 
SIGNAL  ON  MOOULATliO  MICROMAVE  BRIDGE 


FIG  8-2 


i 

I  AFTER  FREQUENCY  SHrT  Af 


I  VECTOR  REPRESENTATION  OF  SIGNALS 
I  BRIDGE 


^  VSWR 


CURVES  USED 


PIG  8-4 


l»0SIT10»i  OP  MN,CM 


NOTES:-  j 

0  OiSTANCE  OS  ALONQ  CURSOR  QVE8  VSWR  CORRfiSRONOiNQ  TOAf  | 

0  DiSTANCE  SA  AROUND  OUTSIDE  CIRCLE  GIVES  ^4  OORRtSRONONG  TD  ] 

®  POSITION  OP  PONT  C  IS  SET  BY  VSWR  AT  RUONANCi  | 

I 

AND  GOVERNS  SIZE  OP  t^EDANCE  CIRCLE  I 

i 

t 

0  NTERSECnON  OP  IMrsOANCE  CIRCLE  AND  CONSTANT  X  ORCLt  GIVU  ! 

POSmON  OP  POINT  S .  AND  CONSEQUENTLY  POSmCN  OP  CURSOR  I 


USE  OF  SMITH  CHART  IN  ANALYSIS  OF  BRIDGE 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  CURVES  FOR 
BRIDGE  RESPONSE 


EFFECT  OF  INACCURACY  IN  SHORT-CIRCUIT  SETTING  ON 
BRIDGE  CHARACTERISTIC 
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(a)  RANGE  A-0 


TRACES  OF  DECAY  OF  IONIZATION  USING  DIFFERENT  RANGES  OF 
CALIBRATION  CURVE  (SEE  FIG  8-9) 
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ELECTRON  CONCENTRATION  No,  CM" 
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TRACE  OF  DECAY  OF  IONIZATION  IN  ARGON 
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DECAY  OF  KDNIZATION  IN  CESIUM  VAPOUR 
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DECAY  OF  ELECTRON  CONCENTRATION  IN  CESIUM  VAPOUR 
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G  8  17  PLASMA  TUBE  CAPSULE  BREAKER  AND  CESIUM  BATH 
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CHAPTER  9 


ATOMIC  BSAJ^  EXPERIMENT 


by 

A. B. Jolly,  T.Bell  and  D. Balfour 


9.1  GENERAL 

The  objective  of  the  atomic  beam  experiment  is  to  measure  the  total  collisicui 
crojs  section  of  cesium  atoms  with  helium  and  to  extend  these  measurements  into 
the  high  temperature  region  (1500°K  to  2000°K)  relevant  to  the  MPD  power  generation 
experiment.  Results  at  low  temperatures  have  been  described  in  earlier  reports. 

To  make  measurements  at  nigher  temperatures  a  Mk  II  apparatus  was  designed  which 
incorporates  high  temperature  heaters  both  on  the  atomic  beam  source  emd  gas 
iscattering  chambers,  has  increased  resolution  due  to  its  increased  size,  and  is 
capable  of  being  pumped  down  to  lower  base  pressures.  This  apparatus  is  now 
complete  and  during  the  last  fcur  months  the  initial  vacuum  troubles  and  cold  spots 
have  been  overcome  so  that  stable  beams  of  atomic  cesium  are  now  produced  regularl;', 
with  an  automatic  temperature  control  of  t  0.5°C  on  the  source.  By  varying  the 
temperature  of  the  nozzle  up  to  1000°K  the  mean  velocity  of  the  atoms  effusing  from 
the  rozzle  may  be  varied  within  the  range  200  to  400  n/sec. 

To  measure  the  intensity  of  the  cesium,  atomic  beam  hot  wire  detectors  relying 
on  siirface  ionization  of  the  cesium  atoms  have  been  used.  Platinum,  niobium  €md 
tungsten  wires  have  been  tested  and  different  diameters  of  wire  have  been  evaluated. 
In  general  the  platinum  and  niobium  detectors,  owing  to  the  greater  purity  of  these 
materials,  give  lower  background  current  than  tungsten  and  wej?e  used  in  the 
subsequent  experimental  programme. 

Electron  bombardment  was  considered  as  a  possible  means  of  heating  the  source 
to  temperatures  in  excess  of  1000°K  and  tests  were  carried  out  but  persistent 
breakdown  of  the  high  voltage  lead*  at  high  temperature  and  irreg\ilar  heating  of 
the  source  chamber  have  impeded  this  development.  The  alternative  of  using  r.f. 
heating  is  being  considered  but  meanwhile  temperatures  up  to  1000°K  are  being 
achieved  using  only  a  'Pyrotenax*  heating  coil. 

The  Mk  II  apparatus  is  equipped  with  a  moveable  rack  which  can  upport  two 
detector  filaments  and  using  this  facility  a  sm'vey  of  beam  profiles  has  been 
carried  out  enabling  alignment  and  detactor  size  to  be  optimized.  The  experience 
obtained  in  this  way  was  then  used  in  experiments  tc  measure  the  total  collision 
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cross  section  of  cesiun  in  helium,  A  conclusion  from  these  experiments  is  that 
the  accuracy  of  results  is  limited  by  the  base  '.ressurc  obtainable  in  the 
scattering  chamber,  Outgassing  from  the  scattering  chamber  may  be  particularly 
disturbing  when  the  chamber  is  heated  to  hign  temperature  and  present  work  is 
directed  towards  eliminating  this  source  of  error. 

The  initial  measurement  of  the  variation  of  cross  section  with  i.he  velocity 

•  2 

of  the  atomic  beam  produced  the  v  t  trend  that  was  an'icipeted  for  cesium  atoms 
in  air  at  5  ^  10"*  tor*r.  However,  on  increasing  the  pressure  of  the  scattering 
gas,  the  results  were  not  reproducible  and  a  similar  measurement  in  helium 
suggested  that  the  cross  section  was  constant  within  the  velocity  range  200  to 
400  av^sec.  Those  problems  are  still  under  investigation  and  must  be  understood 
before  extending  the  range  of  Cs-He  cross  section  measurements  t'^  higher 
temperatures , 

9.2  KXPERIMSMAL  EgUIPlffiOT 

9,2,1  G-eneral  arrangement 

The  Mk  II  atomic  beam  experiment  is  basically  the  same  as  the  Mk  I  experiment 
which  was  described  in  the  eai'lier  Technical  Summary  Reports  (ikii  63“1  and  64-36) 
and  consists  of  a  vacuum  chamber  containing  an  atomic  source  with  a  slit  from 
which  the  alkali  atoms  effuse  as  a  beam,  a  scattering  chamber  with  entrance  and 
exit  slits  aligned  so  that  the  a tomic  beam  can  pass  through  the  scattering  gas  and 
a  hot-wire  detector  to  measure  the  strength  of  the  atomic  beam  after  it  has  passed 
through  the  scattering  chamber.  The  Mk  II  apparatus  is  an  improvement  on  the 
Mk  I:  it  has  improved  resolving  power  due  to  the  greater  separation  between  source 
and  detector,  is  capable  of  being  pumped  to  lower  base  pressures,  and  can  accommodate 
high  temperature  heaters  both  on  the  exit  noezle  of  the  source  and  on  the 
scattering  chaaiber,  thereby  enabling  the  cross-section  measurements  to  be  extended 
to  much  higher  temperatures.  Fig*  9.1  shows  a  general  view  of  the  apparatus  with 
the  source  in  the  foreground,  the  gas  bottles  (above)  connected  to  the  scattering 
chamber  and  the  differential  screw  for  the  detector  mount  at  the  far  end. 

The  general  arrangement  of  the  atomic  beam  oven  and  superheater  is  shown  in 
Pig,  9,2,  A  glass  ampoule  containing  the  alkali  metal  is  housed  in  the  flexible 
bellows  and  can  be  broken  by  compression  of  the  bellows,  allowing  the  alkali  me+al 
to  drop  down  into  the  resei^voir.  The  reservoir  may  be  isolated  from  the  main 
enclosure  and  from  the  nozzle  by  means  of  a  stainless  steel  valve.  This 
facilitates  removal  .■'f  th3  reservoir  for  cleaning  and  rech-arging  without 
disturbing  the  vacuum  condition  existir,g  in  the  main  enclosure.  The  reservoir, 
ampoule  breaker  and  valve  are  attached  to  the  end  flange  of  the  main  enclosure. 
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th«  flange  seal  being  effected  by  compression  of  either  a  silicone  rubber  *0’ 
ring  or  fine  gold  wire. 

The  alkali  vapoi  r  pressure  is  maintained  above  the  liquid  metal  by 
controlling  the  temperature  of  the  stainless  steel  reservoir  which  has  large 
thermal  capacity  to  counteract  the  smell  variations  in  temperature  due  to 
fluctuations  of  the  input  current.  The  vapour  generated  in  the  r-^servoir  then 
passes  to  the  nozzle  through  a  1  in.  bore  stainless  steel  bakeable  valve. 

9*2.2  Source 

A  rectangular  slit  is  cut  vertically  in  the  end  of  tiie  nozzle  and  has 
dimensions  of  0.015  cm  x  0,5  cm.  Care  was  taken  during  the  welding  operation 
to  ensure  that  Ihe  nozzle  slit  was  not  displaced  from  the  axis  of  the  main 
enclosure, 

A  water  cooled  tube  encircling  the  nozzle  absorbs  practically  all  of  the 
radiated  heat,  thereby  ensuring  that  the  enclosure  walls  are  kept  at  a  reasonable 
temperature.  It  also  prevents  oxides  forming  on  the  enclosure  walls. 

The  mean  energy  or  velocity  of  the  atoms  effusing  through  the  slit  in  tha 
nozzle  is  controlled  by  the  temperature  of  the  nozzle*  the  beam  intensity  is 
deterc-ined  by  the  vapour  pressure  in  the  source  chamber,  i.e,  by  the  temperature 
of  the  reservoir. 

The  number  of  molecules  which  emerge  from  the  slit  travelling  in  a  solid 
angle  d«  at  an  angle  6  relative  to  the  normal  of  the  plane  at  the  slit  is 

dQ  =  ^  ET  cosS  A^  ..,,,  (9*1) 

where  n  is  number  density  of  molecules,  v  is  the  mean  molecular  velocity  la  the 

source  and  A  is  the  area  of  the  slit, 
s 

9.2,3  Scattering  chamber 

For  measurements  of  scattering  cross-section  the  scattering  chamber  (Pig.  9,3)  is 
is  aligned  with  the  nozzle  and  detector,  the  slits  of  the  scattering  chamber 
defining  the  beam  width.  The  pressure  in  the  chamber  is  maintained  by  adjusting 
the  needle  valve  of  the  gas  supply  to  balance  the  rate  of  loss  through  the  slits. 

The  scattering  chamber  consists  of  a  stepped  section  (thermal  choke)  tantalum 
cylinder, having  one  end  closed  and  the  other  erd  wolded  to  a  stainless  steel  tuba 
of  1.5  in.  OD  (0,3?k  in.  diameter).  Holes  are  drilled  through  the  tantalum 
cylinder  near  the  closed  end  to  locate  the  tantalum  slit  mountings  which  are  used 
for  adjustment  of  tha  slit  width. 
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A  centrally  mounted  0.25  in.  boro  thin-walled  stainless  steel  tube  is  used 
for  sampling  the  gas  pressure,  which  Is  measured  by  an  ionization  gauge.  Vaa'xous 
spectroscopically  pure  inert  gases  may  be  introduced  to  the  system  from  glass 
h./ttle.«  attached  to  a  5  litre  gas  reservoir.  The  gas  reservoir  is  thoroughly 
outgassed  to  attain  a  pressure  of  10”®  torr  before  Introduction  of  the  scattering 
gas.  The  design  also  provides  for  the  Introduction  of  alkali  metals  so  that 
scattering  cross  sections  can  be  obtained  for  alkali  vapours  and  inert  gas  mixtures. 


9 .2 .4  Surface  ionization  detector 


A  full  desc"iption  of  the  surface  ionization  detector  was  given  in  an  earlier 
report  (IRD  64-36),  but  for  convenience,  the  diagram  is  again  shown  in  Fig.  9.4. 

The  principal  parts  of  the  detector  are  a  thin  filament  heated  to  approximately 
1 000°K,  ensircied  by  an  electrode  having  a  slot  cut  in  line  with  the  filament  for 
admission  of  the  beam.  The  electiode  is  held  at  a  potential  of  -50  volts  with 
I'espect  to  the  fileunent.  For  a  detector  filament  presenting  an  area  to  the 
beam  at  a  distance  r  from  the  slit,  the  number  of  atoms  intercepted  per  second  is 


dQ 


det  4»'r 


A  A^ 
s  d 
7~ 


n  V 


pA.  A 

=  1.12  X  10^^  r^Viff  .  (9*2) 

where  p  is  the  vapour  pressure  in  the  source  expressed  in  torr,  M  is  the  molecular 
weight  in  grams  of  the  beam  particles  and  T  is  the  temperature  of  the  nozzle. 

If  every  particle  incident  on  the  detector  is  registered  as  a  unit  charge 
then  the  maximum  current  attainable  can  be  calculated.  Thus  for  cesium  at  364°K 
and  a  corresponding  vapour  pressure  of  3  x  10"*  torr  the  maximum  current  acceptable 
in  the  Mk  II  apparatus  by  a  niobium  wire  of  0.010  in,  diameter  is  only  2  x  10"^^ amp. 

Neutral  atoms  colliding  with  the  filament  will  normally  be  reflected,  but 
xmder  certain  conditions  may  be  re-evaporated  as  positive  ions  after  a  statistical 
time  interval.  This  so  called  'residence  time*  is  a  function  of  the  surface 
toaperatore  and  work  function  of  the  filamesit  material,  and  the  ionization  potential 
I  of  the  incident  atoms.  The  ratio  of  the  number  of  re-evaporated  ions  n^  to  the 
number  of  re-evaporated  neutral  atoms  n  is  given  by 

r  =  . 

where  k  is  Boltzmanns  constant  (8.5  ^  10"*  eV/°K)  axid  T  is  the  surface  temperature 
in  °K,  It  is  therefore  desirable  that  ((^  -  l)  >>  kT  to  ensure  that  virtually  all 
neutral  atoms  colliding  with  the  detector  filament  are  re-evaporated  and  collected 


as  an  ion  current. 


In  general  it  is  necessary  to  heat  the  filaaent  to  ■''i000°K  (--0.1  eV)  for 
re-emiss j on  to  take  place,  so  that  the  work  function  of  the  filament  material 
should  exceed  the  ionization  potential  of  the  incident  atoms  by  0,5  aV  for 
efficient  surface  ionization  detection.  Cesium  has  an  ionization  potential  of 
3.9  so  that  any  metal  having  a  work  function  of  >  U.U  eV  would  be  a  possible 
filament  material  f or  a  surface  ionization  detector.  The  majority  of  such  metals 
contain  traces  of  alkali  metals  (usually  sodium  and  potassium)  which  give 
troublesome  background  signals;  niobium  and  platinum,  however,  give  low  background 
and  noise  levels  even  when  heated  to  iidOO  K,  Initially  a  tungsten  filament  was 
used  but  with  limited  success  because  of  the  presence  of  thorium  oxides  leading  to 
high  thenaionic  emission  at  the  working  temperature  of  approximately  1200°K.  For 
the  detection  of  sodium  (I  =  5.13  eV) ,  it  is  necessary  to  use  oxygenized  tungsten 
wire  whose  effective  work  function  is  raised  to  about  6  eV. 

Two  detectors  are  provided  in  the  equipment  to  compare  varioi s  filament 
materials  and  to  test  the  resolution  of  the  apparatus.  These  are  mcimted  L.ijacent 
to  each  other  on  a  vertical  plate  which  may  be  accurately  moved  transversely  across 
the  axis  of  symmetry  of  the  apparatus  by  means  of  a  differential  screw,  sliding 
vacuum  seals  and  a  clock  gauge  ( 0,001  in./division). 

Typical  ion  currents  are  10"*^  A;  because  of  the  high  input  impedance  of  the 
current  amplifier  (i0^®n)  the  collector  insulation  must  exceed  10^ ^0  if  a  potential 
of  1  mV  is  to  be  developed  across  a  resistance  of  iOO  MO,  The  use  of  a  Vibron 
vibrating  capacitor  modulator  in  the  ion  current  amplifier  gives  excellent  aero 
stability  (not  more  than  0,1  mV  drift  in  12  hr)  coupled  with  high  input  impedance. 
Simplified  diagrams  of  the  collector  and  filament  circuits  are  shown  in  Figs,  9.5 
and  9.6  . 

9,3  EXPERIMENTAL  RESULTS 

9»3.i  Beam  profiles 

The  movable  detector  mount  in  the  Mk  II  apparatut'  has  made  it  possible  to 
study  beam  prefile  and  the  provision  of  two  separate  detectors  enables  different 
filament  materials  and  wire  diameters  to  be  compared.  Typical  beam  profiles  €ure 
shown  in  Fig.  9.7  where  different  materials  and  different  diameters  of  filament 
were  used  in  the  detector.  The  intensity  scale  has  been  normalized  to  illustrate 
the  effect  of  resolution  on  the  beam  profile;  the  slight  asymmetry  is  due  to 
imperd'ect  alignment  of  the  collimator  slits.  The  ideal  profile  is  shown  for  a 
detector  of  infinite  resolving  power.  As  the  diameter  of  the  detector  filament 
is  decreased  and  hence  the  resolution  increasea  the  profile  tends  towards  the  ideal 
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shape. 

Pig,  9,8  shows  profiles  for  an  unattenuated  beam  and  for  a  beam  that  has 
been  attenuated  by  gas  in  the  scattering  ctiamber.  Collisions  with  gas  molecules 
result  in  broadening  and  attenuation  of  t^e  scattered  intensity  profile.  Whan 
the  scattering  chamber  is  evacuated  the  degree  of  broadening  of  the  beam  profile 
Ab  is  a  measure  of  the  scattering  occurring  with  the  residual  gas  within  the 
vacuum  enclosure.  For  cross  section  determinations  the  value  of  Lh  at  half 
maximum  intensity  was  not  allowed  to  exceed  two  filament  diameters  and  the 
detector  position  was  adjusted  to  record  the  maximum  intensity. 

9, 3*2  Cross  section  measurements 

The  method  of  measuring  total  collision  cross  sections  using  the  present 
experimental  eorangement  has  been  described  in  earlier  reports.  Essentially,  the 
method  consists  of  measuring  the  unattenuated  beam  intensity  and  then  the  attenuated 
intensity  of  the  beam  after  it  has  passed  through  the  scattering  gas  at  known 
pressure  and  temperature.  If  and  I  are  the  intensities  of  the  unattenuated  and 
attenuated  beams  respectively  then 


.  (9.4) 


where  the  brackets  imply  a  mean  value,  1  is  the  effective  scattering  lengtlj  and  X 
is  the  mean  free  path  for  e  beam  a  tom  in  the  scattering  gas.  Then  as  shown  by 
Ramsey'*  : 
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where  a  is  the  mean  effective  collision  cross  section,  p  is  the  gas  pressure  in 
dynes/cm^ ,  k  is  the  Boltzmann  constant,  is  the  gas  temperature  in  °K  and  the 
values  of  g(Z)  where  Z  =  have  been  tabulated  by  Rosin  and  Rabi^  (T,  is  the 

beam  temperature,  amd  are  the  masses  of  the  beam  and  gas  atoms  respectively). 


In  practice  beam  attenuation  is  measured  for  different  values  of  scattering 
gas  pressure  and  the  cross  section  is  then  deduced  from  the  slope  of  log  I  versus 
p  plots.  Pig,  9,9  shows  a  plot  of  log  I  versus  p  for  a  cesium  beam  scattered  by 
xenon  for  a  beam  temperature  of  50o'’k.  In  this  experiment  a  niobium  wire  detector 
of  0,010  in,  diameter  was  used  with  a  resolution  of  ~  1  minute  of  arc.  The  value 
of  collision  cross  section  deduced  from  the  slope  of  this  graph  is  2850  x  10"*®  cm^  . 
Later  measurements  of  this  cross  section  gave  lower  values  similar  to  that  obtained 
in  the  previous  apparatus  but  foi'  these  the  gas  was  less  pure.  The  result  for 
cesium-xenon  is  likely  to  be  realistic  since  the  errors  due  to  low  angle  scattering 
are  considerably  less  than  for  the  case  of  cesium-helium  (described  below). 


9,6 


A  alniilar  tceasurement  of  the  cesium-helium  scattering  cross  section  resulted 
in  the  linear  plot  of  log  I  p  shown  in  Fig.  9.10.  The  value  of  collision  cross 
section  deduced  from  these  results  was  y(  10“‘®  cm^  at  a  resolution  of  0,9 
minutes  of  ‘c.  This  compares  with  total  collision  cross  section  of  1 62  x  10“^®cm^ 
measured  by  Rosin  and  Rabi*  using  a  resolution  of  1  minute  of  arc.  However, 
previous  measurements  in  the  Mk  I  apparatus  and  also  by  Forior^at  higher  resolving 
powers  yielded  much  larger  values.  These  results  imply  that  the  resolving  power 
must  be  increased  to  yield  a  realistic  value  of  total  collision  cross  section. 

It  is  therefore  intended  to  repeat  the  expieriments  using  detector  wires  of  smaller 
diameter. 


9.3.3  Velocity  dependence  of  collision  cross  section 

If  the  nozzle  temperature  be  increased  from  150°C  to  750^C  the  beam  atoms 
effuse  out  of  the  slit  with  an  energy  distribution  determined  by  the  temperature 
of  the  nozzle  and  the  mean  velocity  may  be  varied  within  the  range  200  to  400  nv^sec. 
Since  in  general  the  collision  cross  section  decreases  as  the  energy  of  the  incident 
particles  increases  the  total  number  of  particles  reaching  the  detector  increases 
with  beam  temperature.  Fig,  9.11  shows  the  variation  of  beam  intensity  at  the  detector 
detector  as  a  function  of  beam  temperature ,  In  this  experiment  air  at  a  pressure 
of  5  X  10"*  torn  was  the  scattering  gas  and  a  niobium  detector  ( 0.010  in,  diameter) 

%  18  used.  From  this  curve  it  is  possible  to  deduce  the  velocity  dependence  of  the 
cross  section  according  to  the  argument  which  follows. 


The  attenuation  of  a  well  defined  molecular  beam  in  passing  through 

a  scattering  path  (l)  is  related  tc  the  effective  mean  free  path  of  the  average 
beam  molecule  (velocity  V)  in  a  Maxwellian  gas  by 
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where  a  is  the  total  effective  cross  section. 
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Put 


<7  =  k,  V 
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where  x  denotes  the  velocity  dependence  of  cross-section, 
potentials 

X  =  k2 
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exp 
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Assuming  attractive 

.....  (9.9) 
.....  (9.10) 


so  that 
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(9.11) 


If  the  value  of  x  ia  unlcnoim  then  a  plot  of  log  (log  l)  ~  log  V  will  give  a 

straight  line  with  a  slope  x.  The  value  of  log  is  constant  and  so  may  he 

disregarded.  Hence  the  velocity  dependence  of  cross-section  may  be  directly 

detenained  from  e  xperiasent,  if  the  scattering  process  is  entirely  influenced  by 

the  long  range  Van  der  Waal  attractive  potentials  When  the  results  of  Fig,  9.11 

were  plotted  in  the  above  manner  (Fig.  9.12)  they  showed  an  inverse  0.37  pocer 

i 

velocity  dependence.  This  is  consistent  with  the  dependence  deduced  from 
the  attractive  potential  term  in  the  Lennard-Jones  model. 

Subsequent  measui'ement"  in  helium  have  not  shown  such  a  strong  velocity 
dependence  in  this  temperature  range.  However,  these  are  preliminary  results  and 
further  experiments  are  being  undertaken  to  clarify  this  problem, 
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MICROWAVE  STUDISS  OF  THERMALLY -EXCITED  PUSMAS 


by 

J.  A.  Haines 


10.1  THERMAL  IONIZATION  EXPERIMKITT 


The  object  of  the  thermal  ionisation  experiment  is  to  produce  helium-cesium 
mixtures  in  a  heated  tantalum  tube  at  temperatures  in  the  range  2500°  to  2000°K 
and  to  measure  by  microwave  transmissior.  the  resulting  electron  concentration  (and 
hence  degree  of  ionization)  and  electron  collision  fi*equency.  The  experiment  is 
also  regarded  as  a  laboratory  facility  on  which  to  test  the  microwave  diagnostic 
technique  which  will  eventually  be  used  for  measurements  on  the  MPD  closed-loop. 


In  the  previous  Technical  Summary  Report  (IRD  64-3b)  experiment  was 


described  in  which  the  cesium-helium  mixture  was  heated  in  a  tantalum  tube  to 
2000°K  and  the  shift  in  phase  of  microwaves  propagated  through  the  tube  was  measured 
as  a  function  of  temperature  in  the  range  1500°  to  2000°K.  From  those  measurements 
the  average  electron  concentration  was  deduced  using  the  relationship 


n  =  "*  bdi  (x—  y  ) 


(10.1) 


where  u  is  the  radian  frequency  of  the  microwaves,  64>  is  the  phase  shift  in  radians, 
m  the  mass  of  an  electron,  e  the  electronic  charge  and  c  the  velocity  of  light. 

The  maximum  electron  concentration  detected  was  approximately  5  x  10^°/cm^:  this 
level  of  ionization  is  a  factor  of  100  less  than  the  value  calculated  from  the 
Saha  Equation  on  the  assumption  that  the  gas  pressure  in  the  tantalum  tube  is 
determined  by  the  vapour  density  in  the  cesium  bath.  However,  in  this  experiment 
the  apparatus  was  first  filled  with  helium  and  the  cesium  bath  was  then  heated  to 
raise  the  cesium  vapour  pressure.  In  practice,  diffusion  of  cesium  through  the 
helium  and  through  the  connection  between  bath  and  chamber  is  so  slow  Uiat  the  time 
for  the  cesium  pressure  to  become  constant  throughout  the  apparatus  is  very  large; 
this  condition  may  not  have  been  fulfilled  in  the  experiment.  In  the  next 
exp^eriments  the  cesium  vapour  pressure  will  be  established  before  helium  is 
introduced. 

On  opening  up  the  apparatus  after  the  cesium-helium  experiment,  it  was  found 
that  the  walls  of  the  container  were  coated  with  an  oily  black  substance  and  it 
was  not  again  possible  to  obtain  a  good  vacuum  even  on  baking  to  400°C,  until  the 
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MODIFIED  ASSEMBLY  TO  SUPPORT  THE  TANTALUM  TUBE 


apparatus  had  been  stripped  coapletely  and  scrubbed  with  water  and  determent. 

It  is  believed  that  this  black  coati’*  is  cesius  hytiroxide  formed  when  the  cesium 
coating  on  the  walls  is  exposed  to  moist  air. 

A  more  serious  consequence  of  this  experiment  was  that  the  0,25  iri*  diameter 
tantalum  tube  (0,010  in.  wall  thickness)  fractured  on  cooling.  The  apparat  jis  had 
been  tested  using  a  more  robust  1 .25  in.  diameter  tantalum  tube  which  pennitted  two 
modes  of  microwave  propagation.  To  obtain  interpretable  results  in  which  only 
one  mode  was  dominant  it  was  necessary  to  use  a  tube  of  less  than  0,25  in.  diameter. 

At  high  temperaturn  this  smaller  tube  was  strained  through  inadequate  provision  for 
differential  expansion  forces.  Originally  the  tube  conductor  was  supported  rigidly  on 
on  tantalum  supports  attached  to  a  copper  conductor  rod,  with  the  differential 
expansion  between  the  tube  and  rod  supposedly  taken  by  a  bellows  joint  at  the  end 
of  the  rod.  In  practice,  the  teuitalum  tube  (embrittled  by  temperature  cycling  and 
very  weak  at  2000^K)  was  incapable  of  supporting  the  stress  developed  across  the 
bellows  and  fractured  in  consequence. 

1 0 .2  NEW  DESIGN  OF  SQUIPMENT 

To  eliminate  these  troubles  the  apparatus  has  been  completely  redesigned. 

The  bellows  have  been  removed  and  flexible  copper  leads  carry  the  heating  current 
to  one  end  of  the  tantalum  tube,  the  opposite  and  being  clamped  to  the  end  flange. 

Pig.  10,1  shows  the  modified  structure.  At  the  free  end  of  the  tantalum  tube  two 
stainless  steel  pegs  have  been  fixed  to  the  tantalum  blocks  with  insulating  ceramic 
sleeves  fitted  as  shown.  These  sleeves  rest  on  stainless  steel  support  rods  so 
that  they  can  move  laterally  to  take  up  the  differential  expansion  between  the 
supporting  structure  and  the  tantalum  tube. 

The  assembly  was  tested  by  placing  it  in  a  glove  box  and  heating  it  to  1800°K 
in  £Ln  atmosphere  of  helium  to  simulate  the  experimental  working  conditions.  The 
tube  expanded  laterally  without  any  twi3tir.g  or  bowing  and  the  complete  apparatus 
is  again  being  assembled  in  preparation  for  further  thermal  ionization  experiments. 

There  was  some  difficulty  in  obtaining  a  replacement  tantalum  tube  of  lees 
than  0*25  in,  internal  diameter  and  with  0.020  in.  wall  thickness.  For  a  9  in, 
tube  Accles  and  Pollock  quoted  15  weeks  delivery  time  and  a  cost  of  £150.  The 
Tungsten  Manufacturing  Company  produced  two  seam-welded  tubes  fabricated  from  0,020  in, 
sheet  tantalum  but  these  did  not  have  a  circular  section,  Murex  Ltd  finally 
pro^'uced  a  suitable  tuba  by  drilling  out  the  bore  as  specified  and  then  machining 
away  the  walls  to  0,020  in,  diameter  (Fig.  10.2).  This  particular  tube  has  had 
the  ends  threaded  so  that  it  can  be  screwed  into  the  end  blocks.  In  this  new 
design  the  microwave  horns  will  be  welded  directly  to  the  tantalum  blocks. 


UPPER  SEAM  WELDED  BY  ELECTRON  BEAM  AT  IRD.MDDLE  SEAM 
WELDED  BY  TUNGSTEN  MANUFACTURING  CO.  LOWER  DRILLED 
FROM  SOLID  ROD  GROUND  TO  CORRECT  WALL  THICKNESS 
AND  THREADED  BY  MUREX  LTH 


TANTALUM  TUBES  OF  0  020in.  WALL  THICKNESS  AND  O  23  in.BORE 


Firt-inp 
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REAT  THAjNSFSE  AND  PRESSURE  LOSS  CALCLRJiTIONS  FOR  He-He  RECUPERATIVE  REAT -EXCHANGER 


T,  Archbold 

DESIGN  DATA 

Shell  side  He  inlet  teiLperature  1690°K 

He  outlet  temperature  550*^K 

He  mass  flow  5«75  gw/aeo 

He  average  pressure  1  atm 

Tube  side  He  inlet  temperature  298*^K 

He  outlet  temperature  1438°K 

He  mass  flow  5»75  gja/sec 


A2  GEOMETRY 

The  heat  exchanger  is  of  the  crosa-counterflow  multi-start  helical  coil  type, 
The  shell  side  fluid  flow  through  an  annular  apace  containing  the  coiled  tube 
bvmdle.  Pig,  A,1  shows  the  bundle  details  (excluding  bundle  length).  The  coil 
diameters  and  number  of  starts  per  coil  have  been  chosen  to  give  equal  lengths  for 
every  tube,  and  the  same  helix  angle, 

A3  HEAT  TRANSFER 

A3.1  Shell  side 

The  angle  of  the  coiled  tubes  to  the  overall  direction  of  gas  flow  is  given 


tan  6  = 


where  the  symbols  are  explained  in  the  Nom^nclatm-e  at  the  end  of  the  Appendix, 

The  ceil  diameter  D  is  defined  as 

n 


D  =  2n.F^ 
n  t 

w,2nP, 


so  that  tan  9  - 


=  2v 


6.28 


Hence  S  =  81° 


Thus  it  is  Justifiable  to  consider  the  shell  side  flow  as  normal  cross-flow, 
is  also  reasonable  to  treat  the  coiled  bundle  as  if  it  were  a  banlc  of  straight 
the  curvature  being  perpendicular  to  the  line  of  flow. 

Because  adjacent  tube  layers  are  wound  in  opposite  sersas,  all  forms  of 
in-line  tubs  arrays  with  exis.  in  the  bundle.  The  probable  effect  of 

will  be  to  generate  heat  transfer  rates  greater  than  that  f or  a  square  in-line 
array,  however,  the  buiidle  will  be  designed  using  th.  available  data  for  an 
equivalent  square  in-line  euray  (Fig.  A.2). 


Mininum  flow  area 


^min  expressed  by  Smith^ 


as  follciTs: 


A  . 

min 


V 

4t 


(V-V) 


+  1  + 


p 


where  t  =  F  and  P  =  P 

t  X 

D  =  11 ,1  in.  =  28.2  cm 
o 

=  4.26  in.  =  10.82  cm 

d  =  0*682  in.=  1.73  c® 

t  =  p  =  2.16  cw 

Hence 

=  '27.4 

BquivaT  ent  diameter  i.s  defined  as 

^  _  4  (min.  flow  area)  N.Pj^ 

e  ———————  — — - 

vHeat  transferring  area; 

4  X  127.4  X  2.16  N 

""  1330  *  N 

=  0.828  cm 

5  D 

(T  =  847°C) 

3.73  X  0.828  X  10^ 

^  “  127.4  X  0.0453 

=  77.3 
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Using  the  data  of  Pig.  A. 2  it  ia  clear  that  the  known  Re  range  must  be 
considerably  extrapolated  to  reach  the  necessary  values  of  Re.  As  explained  in 
Kays  and  London’  these  data  were  obtained  using  a  transient  technique  that  is 
reliable  only  in  the  turbulent  region.  This  means  that  the  transient  range, 
indicated  by  a  broken  "^urve,  is  not  accurate.  Consequently,  the  chain  dotted, 
extrapolated  cui've  can  only  serve  as  an  indication  of  the  heat  transfer  data. 
Nevertheless,  it  is  considered  that  a  pessimistic  value  will  be  so  obtained, 
probably  resulting  in  an  overdesigned  tube  bundle. 


The  poor  quality  of  this  data  excludes  the  possibility  of  calculating  the 
tube  bundle  as  a  sequence  of  short  sections  with  different  temperatures. 

From  Fig.  A, 2,  when  Re  =  77*3 

2j 

St.  Pr  ^  3  =  0.031 

at  T  =  847°C,  Pr  =  0.74 

St  =  0.0385 
hA 


Nn  St  = 


min 


G  Cp 


•  K  -  0«0385  X  5.75  x  1.242 
8  127-4 

=  2.l6  X  10’^  oal/sec  om^  °C 

A3. 2  Tube  side 

The  tube  bundle  contains  18  tubes,  so  that  eaoh  carries  a  helium  flow  of 


=  0.31 9  ^0/ sec 

This  will  be  at  an  approximate  pressure  of  1  atm  and  have  a  maximum  temperature 
of  1165°C. 

For  He,pii65%  =  '  gn/um^ 

So  the  maximum  velocity  in  the  pipes  will  bo 

„  _  0.319  X  4  X  1438  _ _  / 

max  1,785  X  10'^  X  X  X  {l.26)*^  x  273  ° 

=  75-3  H/aec 

The  velocity  o;’  sound  in  He  at  1  atm  and  1438°K  is 
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sec 


^  dorn.s/o.’ 


1 .785  X  10'“*  X  2] 

p  =—  u35 - ' 


go/  cm''^ 


_  J_  ri.63  X  14.7  X  10^  X  10'*  X  1438 
100  \  1.45  ’<  1  .785  X  273 

C=2.2xl0^iiV^  sec 

Thus  Mach  No.  in  tubes  =  M  =  _  0.0342  it  is  therefore  possible  to 

max  2.2  x  lO'’  ^ 

neglect  the  effect  of  M. 

Reynolds  No.  for  a  circular  pipe  is 
Re  -  -  ^ 

"  V  Dp 


where 


=  0,319  gw^sec 
=  1  .26  cm 

=  0*042  X  10“^  poise 
_  4  X  0*319  X  10^ 

“  ir  X  1  .26  X  0.042 


Re  =  768 

Thus  the  helium  flow  lies  well  within  the  laminar  region. 

For  the  laminar  flow  of  fluids  in  straight  smooth  pipes,  Hausen^  recommends 


the  relation 


Nu  =  3*65  + 


0.0668 


Pe  D 


1  +  0„045 


>e  dU  J  (‘^w 


For  laminar  flow  in  coiled  pipes,  no  data  a.r-9  known.  It  is  therefore 
proposed  to  use  the  Hausen  form,  modified  by  the  Jeschke  factor  to  include  in  some 
degree,  the  effect  of  the  tube  coiling.  It  should  be  noted  that  the  Jeschke 
factor  has  been  demonstrated  only  for  the  turbulent  region. 


Then 


. 

Nitj,  -  Nu  1  +  3.5  f- 


Provioua  calculations  indicate  that  L  has  a  value  of  about  11  feet,  and  T  —994  K 

'  w 


Hence 


Pe  =  Re.Pr. 

=  0.76  (T^  =  59  fc) 

Re^  =  768 

Re  ^  768  X-  0.76  =  583 


D  =1 ,26  cm 


L  --  11  ft.  =  335 


=  0.042  X  10-^  poise,  =  593°C 
=  0.0443  X  10’ 2  poise,  =  721  °C 


Nu 


=  20.8  cm 


3.65 


^  0.0668  X  583  x  1.26/33^ 
'  ,  .  0.045  ' 


0  *1  4 


3.65  +  0.136 


1 

1.0067 


Ku  =  3.76 


3.76  1 


+ 


3.5  X  1.26 

20.8 


4.58 


Nu  =  h  De/K 


hq. 


4.58  X  6.85  X  10'* 

1 .26 


=  2,49  X  10“^  cal/ssc  cm^  °C 


A3. 3  He&t  transfer  coefficient 

The  overall  heat  transfer  coefficient  U  is  given  by 


where  t  is  the  tube  wall  thickness,  and  is  the  thermal  conductivity  of  Nimonio  75. 
h^  =  2.16  X  10”^  cal/s®c  cm^  °C 
h^  =  2.49  X  10*^  caVsec  cm^  C 
t  =  0.234  cm 


A.5 


Also, 


K  =  0.0622  cal/sec  cai'^C 

w 

i  -  1  1  0.234 

U  “  X  10"'^  *  2.49  X  10'“*'  0.0622 

=  463  +  402  +  5 *76 

~  =  869  aec  cm^  °C/cal 

1  S.  T 
U  "  Q 

S  =  total  heat  transfer  area 
=  1 8  X  jr  X  i  .49  L  ciB^ 

T  =  252°C  (constant  in  recuperator) 

Q  =  Heat  exchanger  power 

=  5*75  ^  1.242  X  1140  cal/sec 

18  X  ff  X  1 .49L  X  252 

5.75  X  1.242  X  1140  "" 

L  =  333  cm 
L  =  10.9  ft 

Considering  the  inner  coll,  the  length  of  each  pipe  is 

L  =  m +  n*p^ 

D  =  12.96  cm 
n 

n  =  3 

p  =  2.16  cm 

m  =  no.  of  turns  per  pipe  in  layer 
“  ^^y^joTx  12.96)2  +  9  X  2.162 

m  =  8.05 

Core  bundle  length  y  =  m.n.p 

=  8.05  X  3  X  1.16  cm 

y  =  52  cm 

Considering  the  highly  in.accurete  data  used  in  this  calculation,  it  is  recommended 
that  the  coil  bundle  should  be  built  at  least  2^  longer. 
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Thus  the  coil  will  be  at  least  65  cn  long. 


A4 


PRESSURE  LOSSES 


A4.1  Shell  side 


From  Kays  and  London^,  the  total  shell  side  pressure  losses,  including  the 
inlet,  outlet  and  case  friction  effects  are: 


AP  = 

s 


min  1 


/  min  1 ^ 

+  ) 

*f 


^  - 1 
P7 


+  4  f 


£1-1 
“  j 

ik  fs 


C  =  5*75  gs/ sec 

A  =  127*4  cm^ 
min 

Pi  =  1.785  X  10~^  X  =  2.88  X  10“^  ffo/ om^ 

=  533  cm* 

pj  =  1.785  X  ICT^  X  =  S.83  X  10"®  gn/cm* 

f  -  0.045  (see  Fig.  A.2) 
y  =  65  cm 
D  =  0.828  cm 

p  =  5.85  X  10"®  gnv/cm* 

NB«  All  densities  eure  calculated  for  1  atm 

AP  =  220  dyne  a/ cm* 

8 

.  .  AP^  =  0.0032  psi 
A4.2  Tube  side 

Pressare  loss  in  a  tube  with  inlet  fluid  temperature  Ti  K  and  outlet  T2K,  with 
nnifonn  heat  flux  is  given  by 


For  streamline  flow,  in  curved  pipes,  the  normal  friction  factor  derived 
for  straight  pipes  must  be  increased  by  a  factor  developed  by  Drew  (in  McAdams^). 

The  average  of  this  factor,  weighted  according  to  the  number  of  pipes  per 
layer  =  1 .84 


.*.  f  =  1 .84  K  16/768 


f  =  1 .84  X  16/Re 
Re^  =  768 

f  =  0.0383 

G  =  0.319  fW sec 

=  1,785  X  10"'*  X  =  5.62  X  10"*  gis/cm^ 

=  1 .245 

L  =  333  cm 

Tj=  1438°K->_ 

T  =  866°K 

T,  =  298‘*K 

AP  =  2.52x10*  dynea/cm^ 

AP^=  0.365  psi 

A5  PHYSICAL  PROPERTIES  OF  Jle 
Cp  =  1 .242  cal/gia®C 

r  1 .63 

p  =  1.785  10"'*  gn/om^  (273°K  and  1  Atmos) 

p  -  McAdams*  (pp  468-469  Table  A-19  and  Fig.  A-6) 

K  -  Eckert®  (up  to  675*^0). 
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NOMSljCMTURS  (APPEND  II ) 


Cp 

D 

e 

"n 

D 

o 

d 

f 

r* 

h 

K 

L 


ts 

U 

n 

Nu 

p 

X 

P 

Pr 

Pe 

Q 

Re 

St 

S 

T 

t 

U 

V 

y 


miniffluffi  flow  area 

shell  side  frictional  eirea 

velocity  of  sound 

spec  heat  at  constant  press 

equivalent  diaiaeter 

helix  diaseter 

diameter  of  layer 

shell  outer  diameter  (inside) 

shell  inner  diameter 

tube  outside  diameter 

fanning  friction  factor 

mass  flow 

heat  transfer  coefficient 

thermal  conductivity 

tube  length 

no.  of  pipes  in  layer 

Mach  niunber 

no.  of  layers 

Nusselt  numbers 

transverse  pitching  of  tubes 

axial  pitching  of  tubes 

pressure 

Prandtl  number 

Peclet  number 

heat  exchanger  power 

Reynolds  number 

Stanton  number 

total  heat  transfer  area 

temperature 

tube  wall  thickness 

overall  heat  transfer  coefficient 

velocity 

bundle  length 


Sreek  symbols 

f  ratio  of  specific  heats 

AP  pressure  differences 

AT  temperature  differences 

p  density 

6  angle  of  tube  to  gas  flow  direction 

p  viscosity 

viscosity  at  weQ.1  temperature 
P^^  viscosity  at  mean  bulk  temp,  of  fluid 


Suffixes 

S  shell  side 

t  tube  side 

wall 

mean  value 


w 


SOURCE  comrrct  exchangers 


HEAT  TRANSFER  AND  FRICTION  FACTOR  DATA 


FIG  A- 2 


EXTKIIML  DISTRIBUTION 


Director,  Advanced  Research  Projects  Agency 
The  Pentagon 

Washington,  D.C. 20301  (2) 

Office  of  Naval  Research 
Power  Branch  (Code  429) 

Washington,  D.C. 2 0360 

Attn:  John  A,  Satkowski  (6) 

Commanding  Officer 

Office  of  Naval  Research  Branch  Office 

Box  39  Navy  ^100  Fleet  Post  Office 

New  York,  New  York  (1 ) 

Commander  John  £,  Raymond 
Office  of  Naval  Research 
Key sign  House 

London  W,1 ,  (1 ) 

Drc  J,H.  Huth 
Department  of  the  Navy 
Bureau  of  Ships, 

Washington  25,  D.C. 

U.S,  Naval  Research  Laboratory 
Washington  25,  D.C. 

Attn:  Technical  Information  Division  (6) 

Wright-Patterson  Air  Force  Base 
Aeronautical  Sratems  Division 
Ohio 

Attn:  Lt«  Robert  Barthelemy  (AFAPL)  (i) 


Defense  Documentation  Center 
Camerson  Street 

Alexandria,  Virginia  22314  (20) 

U.S.  Army  Research  &  Development 
Laboratory 

Fort  Belvoir,  Virginia 

Attn:  Frank  Shields  (SRD-EP)  (l) 

NASA,  Lewis  Research  Center 
21,000  Brookpark  Road 
Cleveland  35,  Ohio 

Attn:  Wolfgang  Moeckel  (l) 

Dr.  B,  Lubarsky  (l) 

Dr.  U.  Kosmahl  (l} 

U.S.  Atomic  Energy  Commission 
Division  of  Reactor  Development 
Direct  Energy  Conversion  Section, 

RD:  AED 

Germantown,  Maryland  (l) 

Attn:  Dr.  R.  Grossman  (B301 )  (1 ) 

Mr.  T.  Brogan 

AVCO  -  Everett  Research  Laboratory 

2385  Revere  Beach  Parkway 

Everett,  Massachusetts  (I) 

Dr.  J.  Cole 

Department  of  Aeronautics 
California  Institute  of  Technology 
Pasadena,  California  (1) 


Air  Force  Office  of  Scientific  Research 
Washington  25,  D.C. 

Attn;  Dr.  Milton  M.  Slawsky  (1) 

U.S.  Naval  Ordnance  Test  Station 
Propulsion  Applied  Research  Group 
China  Lake,  California 
Attn:  Leroy  J.  Krzycki  (Code  4506)  (1) 

Rome  Development  Center 
Rome,  New  York 

Attn;  Mr.  Frank  J.  Mellura  (1 ) 

U.S.  Naval  Ordinance  Laboratory 

NA  Division 

White  Oak,  Maryland 

Attn;  Wallace  Knutsen  fl ) 

Library  {2) 


Mr.  Arthur  Sherman 
General  Electric  -  Valley  Forge 
Valley  Forge  Space  Technical  Center 
Philadelphia  1,  Pennsylvania  (1) 

Dr.  W.D.  Jackson 

ElectriceJ.  Engineering  Department 
Massachusetts  Institute  of  Technology 
Cambridge  39,  Massachusetts  (0 

Dr.  Vernon  H.  Blackman 
MHD  Research  Incorporated 
1535  Monrovia  Street 

Newport  Beach,  California  (1) 


Dr.  Robert  Bustis 

Thermosci  nces  Division 

Stanford  University 

Stanford,  California  (l) 

Dr.  Richeird  Schamberg 

Rand  Corporation 

1700  S.  Main  Streat 

Santa  Monica,  California  (1) 

Dr.  Sam  Naiditch 
Unified  Science  Associates 
826  Arroyo  Parkway 

Pasadena,  California  (1) 


Dr.  W.S.  Emmerich 

Westinghouse  Research  Laboratories 
Beulah  Road,  Churchill  Borough 
Pittsburgh  >5,  Pennsylvania 

Dr.  R.T.  Schneider 
Allison  Division 
CTeneral  Motors  Corporation 
Indianapolis,  Indiana 

Dr.  D.G.  Elliott 

Jet  Propulsion  Laboratory 

Pasadena,  California 


